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ABSTRACT 

Ibis  report  describes  recent  Investigations  of  the  data  obtained  by  the 
Jason  sounding  rockets  during  the  Argus  tests.  These  instrumented  rockets 
detected  electrons  which  were  injected  into  the  geomagnetic  field  by  the 
detonation  of  nuclear  devices  at  high  altitudes.  The  original  analysis  of 
the  data  has  been  published  in  the  open  literature;  the  further  investigations 
reported  here  yielded:  (l)  the  orientations  of  the  detectors  with  respect  to 
a  plane  perpendicular  to  local  magnetic  lines  of  force,  (2)  the  calibration 
of  the  detectors  in  a  planar  flux  of  electrons  as  functions  of  electron  energy 
and  orientation  of  the  detectors  with  respect  to  the  electron  plane,  and  (3) 
a  re-evaluation  of  the  Jason  data  based  on  items  (l)  and  (2)  and  on  recent 
data  on  characteristics  of  the  natural  radiation.  The  conclusions  verified 
the  results  of  the  initial  analysis  regarding  the  location,  width,  and 
stability  of  the  electron  "shell"  produced  by  the  Argus  II  detonation,  as 
well  as  the  presence  of  a  less  intense  electron  flux  in  wide  regions  (wings) 
away  fro*  the  shell.  A  point  of  departure  was  that  the  electron  spectrum  was 
much  more  in  agreement  with  the  fission  beta  spectrum  than  that  resulting  from 
the  initial  analysis.  Within  the  accuracy  of  the  data,  the  spectrum  was  found 
to  be  independent  of  the  spatial  locations  of  the  rockets  or  of  the  passage  of 
time.  The  decay  of  the  flux  appears  to  follow  a  time  dependence  of  •15+.03 
and  a  lower  limit  to  the  injection  efficiency  appears  to  be  0.013. 
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1.  INTRODUCTION 

The  Jason  data  refer  to  the  information  on  the  geomagnetic  trapping  of 
electrons  obtained  by  sounding  rockets  during  the  Argus  tests.  These  in¬ 
strumented  rockets  detected  electrons  which  were  injected  into  the  geomag¬ 
netic  field  by  the  detonation  of  nuclear  devices  at  high  altitudes.  Project 
Jason  was  directed  by  the  Air  Force  Special  Weapons  Center,  and  personnel  of 
the  Physics  Division  of  the  AFSWC  performed  Tine  initial  analysis  of  the 
data.^ 

A  more  thorough  analysis  of  these  data,  however,  was  dictated  by  the 
military  and  geophysical  importance  of  the  Argus  experiment  and  the  restric¬ 
tion  on  the  continuation  of  similar  experimental  investigations.  In  particu¬ 
lar,  analytical  work  was  needed  to  compute  the  body  motions  of  the  rockets  in 
free  flight  in  order  to  correlate  the  counting  rates  of  the  detectors  with 
the  orientations  of  the  detectors.  Also,  experimental  work  was  required  to 
determine  the  response  of  each  detector  as  a  function  of  electron  energy  and 
as  a  function  of  the  orientation  of  the  detector  with  respect  to  the  distri¬ 
bution  of  the  electron  flux.  By  establishing  these  essential  characteristics 

2 

of  the  measurement  system,  and  also  by  utilizing  relatively  recent  information 
on  the  spectral  and  spatial  properties  of  the  natural  radiation,  a  better 
deduction  of  the  results  of  the  Argus  experiment  was  expected. 

This  further  analysis  of  the  Jason  experiment  is  given  in  this  report. 


1.1  REVIEW  OF  EXPERIMENT 

In  the  Jason  program,  nineteen  four-stage  solid-propellant  rockets  were 
launched  from  Wallops  Island,  Virginia,  Cape  Canaveral,  Florida,  and  Ramey 
Air  Force  Base,  Puerto  Rico.  Table  1-1  gives  a  summary  of  the  rocket  program. 


1.  "Project  Jason  Measurement  of  Trapped  Electrons  from  a  Nuclear  Device  by 

Sounding  Rockets,"  by  Lew  Allen,  Jr.,  James  L.  Beavers  II,  William  A. 
Whitaker,  Jasper  A.  Welsh,  Jr.,  and  Roddy  B.  Walton,  Proceedings  of  the 
National  Academy  of  Sciences,  August  1959;  and  Journal  of  Geophysical 
Research,  64,  893,  1959. 

2.  "Energy  Spectrum  and  Angular  Distributions  of  Electrons  Trapped  in  the 

Geomagnetic  Field,"  by  J.  B.  Cladis,  L.  F.  Chase,  W.  L.  Imhof,  and 
D.  J.  Knecht,  Journal  of  Geophysical  Research,  66,  2297,  1961. 

1 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


Launch  coordinates  for  the  Jason  Flights  vere:  Wallops  (75°29'W,  37°50'N)j  Partick  (80°32'V,  28°2 7'R); 

Ramey  (67°7'W,  l8°32'N) 
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Rockets  which  performed  properly  generally  reached  altitudes  of  7 00  to  800 
km.  Six  of  the  rockets  failed  to  operate  as  planned.  Ten  of  the  remaining 
rockets  were  launched  successfully  following  the  second  nuclear  detonation. 

Of  these,  only  the  five  rockets  which  were  launched  from  Wallops  Island, 

Virginia,  traversed  the  electron  shell.  These  rockets  were  launched  at  time 
intervals  after  the  burst  ranging  from  28  minutes  to  about  89  hours. 

The  radiation-detection  instruments  consisted  of  eight  geiger  tubes  with 
various  absorbers  and  collimators.  A  diagram  of  the  instrumentation  package 
is  shown  in  Figure  1-1.  The  detectors  in  Channels  1,  5,  7,  3»  and  6  are 
referred  to  in  the  text  as  the  long  detectors.  Each  of  these  detectors 
utilized  the  same  type  of  geiger  tube,  an  Anton  106-C.  Of  these,  only  the 
detectors  in  Channels  1  and  5  were  rather  highly  collimated,  as  shown  in 
Figure  1-1.  The  detectors  in  Channels  2,  8  and  4  are  referred  to  in  the  text 
as  the  short  detectors.  These  also  utilized  a  single  geiger  tube  type,  but 
it  was  an  end-window  Anton  222R.  All  of  the  222R  detectors  were  highly  colli¬ 
mated,  each  in  the  manner  shown  for  the  detector  of  Channel  2  in  Figure  1-1. 
Roughly,  the  absorbers  were  selected  such  that  the  detectors  in  Channels  7 
and  8  had  thresholds  at  about  440  kevj  the  detectors  in  3  and  4  had  threshold 
energies  of  about  1  Mev,  and  the  detector  in  Channel  6  was  alone  with  a  thres¬ 
hold  energy  of  4.3  Mev.  The  detectors  are  described  more  fully  in  Section  3« 

The  output  pulses  of  the  geiger  tubes  were  sequentially  sampled  for  a 
time  duration  of  about  1/75  sec  with  a  multi-segment  commutator,  and  the  pulses 
were  transmitted  directly  to  ground  stations.  In  the  telemetry  records, 
regularly  spaced  synchronization  pulses  appeared,  followed  by  pedestal-shaped 
signals,  each  of  which  corresponded  to  a  sampling  of  the  output  of  a  particular 
detector.  The  geiger  tube  pulses  appeared  at  the  top  of  the  pedestals.  For 
each  Channel,  the  counts  per  pedestal  were  determined  visually  by  personnel  of 
the  Air  Force  Special  Weapons  Center  and  recorded  as  a  function  of  time. 

These  data,  together  with  the  trajectories  of  the  rockets  and  the  telemetry 
signal  strength  records  received  by  the  six  ground  stations  of  the  Air  Force 
Missile  Test  Center  range  were  used  in  this  study. 
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1.2  GENERAL  METHOD  OF  ANALYSIS 

A  correlation  of  the  counting  rates  of  the  detectors  with  the  orientation 
of  the  longitudinal  axis  of  the  rocket  in  space,  as  described  in  the  following 
section,  revealed  that  the  distribv  tion  of  the  electron  flux  was  symmetrical 
about  a  plane  perpendicular  to  the  local  geomagnetic  lines  of  force.  Moreover, 
the  angular  "spread"  of  the  electron  paths  about  this  plane  (pitch-angle 
distribution)  was  found  to  be  narrower  than  the  angular  apertures  of  the 
collimated  detectors  (half-angle  at  half-height  of  about  15°). 

It  was  unfortunate  that  the  detectors  were  not  collimated  sufficiently  to 
permit  a  measure  of  the  pitch-angle  distributions  of  the  electrons  and  hence 
obtain  the  distribution  of  the  flux  along  magnetic  lines  of  force.  But  an 
initial  experiment  in  an  unknown  area  rarely  yields  all  that  is  desired! 

However,  since  the  distribution  of  the  electron  flux  about  the  normal 
plane  fell  within  the  collimation  apertures  of  the  detectors,  the  total  flux 
and  the  spectrum  of  the  electrons  along  the  paths  of  the  rockets  could  be 
determined  by  relating  the  counting  rates  and  the  response  functions  of  the 
detectors  to  the  orientation  of  the  detectors  with  respect  to  this  normal  plane. 

This  method  of  analysis  yields  a  planar  flux;  that  is,  a  flux  distribution 
such  that  the  electron  velocity  vectors  are  uniformly  distributed  about  every 

point  in  planes  normal  to  geomagnetic  lines  of  force  but  have  no  components 

_2 

out  of  these  planes.  If  a  planar  flux  is  given  in  units  of  number  •  cm 

sec”  ,  the  numerical  value  of  the  flux  equals  the  number  of  electrons/second 

2 

which  traverse  a  spherical  detector  of  1  cm  cross  sectional  area.  Since  such 
a  detector  would  be  traversed  by  the  same  number  of  electrons  per  second  if 
located  in  an  omnidirectional  flux  numerically  equal  to  that  of  the  planar  flux 
and  given  in  the  same  units,  the  planar  flux  evaluated  by  the  analysis  dis¬ 
cussed  below  is  numerically  identical  to  the  omnidirectional  flux.  The  plane 
of  the  flux  (normal  to  geomagnetic  line  of  force)  is  referred  to  in  this 
report  as  the  radiation  plane. 

Now,  if  the  angles  Tj  and  are  used  to  specify  the  orientation  of  a 
detector  with  respect  to  the  radiation  plane,  then  the  corresnonding  geometric 
function  of  the  detector  will  be  defined  by  means  of  the  equation, 
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V  7 


y) 


IfWG^Vj  ,y,  B)dE, 


(1-1) 


vhere 


f(E) 

Gi(7?'T,E) 


corrected  counting  rate  of  gelger  tube  in 
Channel  i,  oriented  at  angles  71  and  /Mf 
to  the  radiation  plane,  ' 

planar  flux  of  electrons  in  units  of  number 
cm”2  ,  sec-1  .  Mev  . 

geometric  function  of  detector  in  Channel  i 
(units  of  cra^  for  detection  of  planar  flux 
of  electrons  of  energy  E  when  oriented  at 
angles  Of  and  77  to  the  radiation  plane. 

Mote  that  G  includes  the  detection  efficiency 
of  the  gelger  tube. 


With  the  counting  rates,  R^ ,  correlated  with  Yj  and  Ojr  as  described 
in  the  following  section  and  the  geometric  function  G(YJ  E)  determined 

experimentally  as  described  in  Section  3,  the  planar  flux  f(E)  was  "unfolded" 
with  a  digital  computer  from  the  integral  equations  (l-l).  The  computer 
program  used  to  deduce  the  flux  is  described  in  Section  4. 

A  digital  computer  was  used  to  correlate  the  counting  rates,  R^,  with  the 
orientation  angles  Yj  and  ^  .  These  orientation  angles  were  computed  for 
each  time  at  which  a  detector  output  was  sampled.  The  read-out  record  from 
the  computer  listed  the  following  information  for  each  sampling  of  a 
detector  output:  the  time  of  the  sampling;  the  counts  per  pedestal;  the 
corrected  counting  rate;  and  the  orientation  angles  Yj  and  ^  .  A  "dead 
time"  correction  of  70  microseconds  was  applied  to  the  counting  rates.  This 
"dead  time”  for  the  detection  system  was  estimated  from  the  time  resolution 
of  the  gelger- tube  pulses  in  the  oscillograms  from  which  the  counts  per  pedes¬ 
tal  were  read. 

The  values  of  R^  used  in  Equations  (l-l)  were  those  which  were  averaged 
over  certain  intervals  of  time  and  angles  Yj  and  .  This  was  done  to 
improve  the  overall  accuracy  of  the  results. 
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2.  BODY  MOTIONS  OF  ROCKET  IN  FREE  FLIGHT 

The  body  motions  of  the  rockets  in  free  flight  were  determined  by  using 
both  the  counting  rates  of  the  detectors  as  a  function  of  time  and  the 
antenna  signal  strength  records.  The  analysis  conducted  to  evaluate  the 
time-dependent  angles  Tj  and  Hj/'  which  gave  the  orientation  of  each 
detector  with  respect  to  the  radiation  plane  is  described  in  the  section 
below.  In  the  following  section,  the  evidence  is  given  that  this  radiation 
plane  was  normal  to  the  local  geomagnetic  lines  of  force. 

2.1  ORIENTATION  OF  DETECTORS  WITH  RESPECT  TO  RADIATION  PLANE 

The  analysis  of  the  motions  of  the  detectors  with  respect  to  the  radi¬ 
ation  plane  will  be  described  with  reference  to  the  diagram  shown  in  Figure 
2-1.  Here,  a  sphere  is  depicted  with  its  center  at  the  center-of-mas3  of 
the  rocket.  The  labelled  dots  on  the  sphere  represent  the  points  at  which 
the  following  vectors  intersect  the  surface  of  the  sphere:  S  is  aligned 
with  the  spin  axis  of  the  vehicle,  A  is  the  constant  angular  momentum 
vector,  B  is  the  local  geomagnetic  field  vector,  C  is  a  vector  in  the  plane 
normal  to  S  and  directed  along  the  center  of  the  angular  apertures  of  the 
detectors  in  Channels  1,  3,  5,  6  or  7  (for  the  detectors  in  Channels  2,  4, 
or  8,  C  is  along  the  projection  of  the  collimator  axis  in  the  plane  normal 
to  S).  As  shown  in  the  diagram,  0  is  the  angle  between  A  and  B,  7 'j  is  the 
angle  between  B  and  S,  and  is  the  angle  between  A  and  S.  Three  lines 

of  nodes  are  also  shown:  n(6,  A)  is  drawn  normal  to  S  and  A,  n(S,S)  is  norial 
to  S  and  S,  and  n(A, S)  is  normal  to  A  and  S. 

The  vector  C  rotates  in  a  clockwise  sense  about  S  with  the  constant 

angular  velocity  6J0( spin  rate).  The  angular  displacement  of  C  will  be 

represented  by  ©  (t)  =  0 vj  t  +  d  ,  measured  from  the  S,  A  plane.  The  spin 
s  s  s 

axis  S  rotates  in  a  cone  about  A,  in  the  clockwise  sense,  and  with  the  constant 
angular  velocity  (precession  rate).  The  angular  displacement  of  the  A,  S 

plane  will  be  designated  by  ©^(t)  =  t  +  d^,  measured  from  the  node 

n(8,  A). 
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From  this  diagram  the  following  relationships  can  be  obtained: 

cos  Yj  -  cos  0  cos  C  -  sin  0  sin  £  sin  ( CD  +  d  ),  (2-1) 

i  P  P 

Of/'  a  t(J  t  +  d  -  where  (2-2) 

f  s  s 

2/(t)  ■  7T  -  tan”^  £  sin  £  ctn  0  +  cos  ^  sin  (l<9pt  +  dp)J 

.  sec  (<Opt  *  dp)  j  (2.3) 

Hence,  the  orientation  angles  /YJ  and  sy'  can  be  determined  as  functions  of 

time  by  evaluating  the  constants  L\)  ,  CO  ,  d  .  d  ,  0  and  ^  . 

p  s  p  s 

Fortunately,  not  all  of  these  constants  axe  mutually  independent.  The  inter¬ 
dependencies  of  these  constants  and  the  methods  used  to  evaluate  them  are 
discussed  below. 


2.1.1  DETERMINATION  OF  PRECESSIONAL  ANGULAR  VELOCITY 


The  constant  CO  was  determined  from  the  telemetry  signal  strength 
P 

information  which  was  recorded  by  a  number  of  ground  stations.  A  rotation 
of  S  in  a  cone  about  A  was  discernible  in  the  record  because  the  signal 
strength  varied  in  accordance  with  the  aspect  of  the  antenna  pattern 
(quadrupole)  to  the  receiving  station.  Hence,  the  period  of  precession, 

Tp  =  2TT/C Ol  ,  was  given  by  the  time  duration  of  a  recurrent  modulation 
pattern  in  the  record. 


In  order  to  avoid  errors  in  this  measurement  due  to  a  rapid  azimuthal 
motion  of  the  line-of-sight  between  the  rocket  and  the  ground  station,  the 
signal  strength  record  of  a  station  which  was  most  nearly  in  the  plane  of 
the  rocket  trajectory  was  used.  The  records  of  six  ground  stations,  located 
south  of  the  rocket  trajectory  and  spanning  a  longitude  range  from  30°  32'  W 
(Patrick  )to  6l°  45'  W  (Antigua),  were  available  for  this  purpose. 


The  examination  of  the  signal  strength  records  revealed  that,  for  every 
flight  except  2027,  "the  precession  period,  Tq,  was  essentially  constant  over 
the  entire  portion  of  the  trajectory  above  the  "sensible"  atmosphere  of  the 
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earth.  The  rocket  on  2027  exhibited  an  appreciable  diminution  of  rotational, 
energy  during  the  course  of  the  flight:  The  precession  period  increased  by 
nearly  1$  and  the  spin  perio1  increased  by  about  20$.  The  re  Is  also  some 
evidence  of  such  a  draft  on  '".light  2024,  but  the  magnitude  of  the  drift  was 
smaller  by  an  order  of  magnitude.  A  rough  calculation  indicates  that  the 
decrease  of  the  rotational  velocity  of  the  vehicle  in  free  flight  was  caused 
by  the  interaction  of  the  geomagnetic  field  with  the  rotating  eonductive 
structures  and  components  of  the  rocket,  dissipating  the  rotational  energy 
in  the  form  of  eddy  currents. 

The  value  of  T  given  by  this  analysis  are  listed  in  the  second  column 
P 

of  Table  2-1. 


TABLE  2-1 


Flight  No. 

Precession 
Period,  Tp 
(sec) 

Spin  Period, 

Ts 

(sec) 

0 

Degrees 

c 

Degree 

2019 

3.793  +  *001 

4.861  +  .oo4 

5  +  2 

89.5 

2021 

30.00  +  .02 

6.654  +  .005 

16.2  +  1.5 

87.1 

2024 

17.41  +  .01 

1.97  +  .02 

20.5  +  1.6 

84.5 

2027 

1.818  to  I.O30 

1.44o  to  1.756 

5  ±  3 

89 

2042 

9.958  +  .003 

4.292  +  .002 

4  +  2 

88.5 

2.1.2.  DETEFWINATION  OF  ROTATIONAL  ANGULAR  VELOCITY 

The  angular  velocity,  of  the  rocket  about  its  longitudinal  axis, 

3,  was  obtained  from  the  "radiation  signal, "  defined  as  the  curve  obtained 

i 

by  plotting  the  counts  per  pedestal  of  a  channel  versus  time.  (Principally, 
Channels  1  and  5  were  used  for  this  purpose.)  By  measuring  N,  the  number  of 
cycles  (including  fractions  of  a  cycle)  of  the  radiation  signal  which 
occurred  in  an  integral  number,  n,  of  precession  periods,  the  value  of  CxJ  „ 
was  deduced  from  the  relation, 

OJ  B  =  (N/2n)6Up.  (2-4) 
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This  equation  was  appropriate  because,  for  every  rocket  flight,  the 
angle  C  was  found  to  be  greater  than  0  as  described  below.  If  had 
been  smaller  than  0,  then  the  equation, 

CJB  =  (N/2n  -  1)  uj  ,  (2-5) 

S  P 

would  have  been  used# 


The  equation^  above  can  be  verified  in  the  following  manner:  Differ¬ 
entiation  of  Equation  (2-2)  with  respect  to  time  give6, 


+  U)v  (  cos  C  -  cos  0  cosTJ  )/sin .  (2-6) 

Since  ^  is  measured  from  the  radiation  plane  while  t  +  dg  is  in  the 
first  or  second  quadrants,  as  shown  in  Figure  2-1,  equals  one-half  the 
angular  frequency  of  the  radiation  signal.  Hence,  the  average  value  of 
over  the  integral  number  n  of  precession  periods  is, 


</f>  =  (1/2)  •  2TT  N/nTp 

=  (l/nTp)  (  (cos^  -cos  0  cos  Yj  )/sin^yJ  dt 

Jo 


By  substituting  in  this  equation  the  value  of  cos  rj  given  by  Equation  (2-l) 
and  performing  the  integration,  equations  (2-4)  and  (2-5)  can  be  obtained. 
(After  substituting  Equation  (2-l)  for  cos  Yj  ,  the  second  term  in  the  inte¬ 
gral  above  can  be  integrated  by  partial  fractions.) 

The  measurements  of  the  spin  periods,  Tp  =  27T /(^j  g,  of  the  rockets  are 
listed  in  the  third  column  of  Table  2-1. 


2.1.3  DETERMINATIONS  OF  0  AND 

In  principle,  three  independent  methods  can  be  used  to  obtain  the  value 
of  4T  )  and  'two  independent  methods  can  be  used  to  evaluate  0.  First,  the 
angle  can  be  evaluated  by  considering  the  vehicle  to  have  the  configura¬ 
tion  of  a  "symmetrical"  top.  Then  if  R  represents  the  ratio  of  the  small  to 
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the  large  principal  moments  of  inertia  of  the  vehicle,  C  is  given  by  the 
equation 

cos  C  =  ft  bjjco  p*  (2-7) 


The  angle  C  can  also  be  determined  by  analyzing  the  signal  strength 
records.  This  method  is  discussed  in  the  following  section.  The  angle  9 
can  be  obtained  in  this  way  by  using  available  geomagnetic  dip  and  declina¬ 
tion  information. 


Finally,  the  angles  and  0,  as  well  as  the  precessional  phase  angle, 
d  ,  can  be  determined  by  measuring  the  time  variation  of  the  period  of  the 

p 

radiation  signal  within  a  precession  period.  Equation  (2-6)  gives  the  time 
dependence  of  one-half  the  angular  frequency  of  the  radiation  signal;  it 
varies  from  a  minimum  value  of 


Tmin  m  -  ^p  Sln  0/sin  (  C  ’  0) 


for  t  +  d  =  -  T^/2, 
P  P  ' 


(2-8) 


to  a  maximum  value  of 

Tm 


max 


=  <•  CV  sin  0/sin  (  C  +  «)> 

b  P 

for  CO  t  +  d  =  ^/2. 

P  P 


(2-9) 


Hence,  a  graph  of  the  time  interval  between  the  peaks  of  the  radiation  signal 
gives  the  precessional  phase  angle,  d^,  and  the  extreme  values  of  0^  • 
Equations  (2-8)  and  (2-9)  can  then  be  solved  simultaneously  to  yield  9  and 

C  - 


•  • 

Actually,  the  extreme  values  of  /jr  could  not  be  measured  sufficiently 
well  to  determine  both  £  and  0  with  good  accuracy.  A  graph  of  some  of  the 
measurements  of  the  period  of  the  radiation  signal  within  a  precession 
period  of  the  rocket  is  shown  in  Figure  2-2.  The  time  marked  0  on  the 
abscissa  corresponds  to  the  time  07**5  seconds  recorded  on  the  data  read-out 
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sheets*  From  this  point  on,  differences  in  the  times  of  occurrence  of  the 
peaks  were  plotted  against  the  time  between  the  peaks  given  by  the  "read¬ 
out”  time  minus  the  sum  of  07^5  nee  and  an  integral  number  of  precession 
periods  (17*41  sec,).  The  scatter  of  the  measurements  shown  in  the  graph 
are  representative  of  the  data  taken  throughout  the  flight. 

The  value  of  £"  given  by  the  signal  strength  analysis  for  flight  2024 
agreed  with  the  value  of  C  given  by  Equation  (2-7)  to  within  1°.  The 
value  R  =  l/90  was  used  for  the  ratio  of  the  moments  of  inertia.  There¬ 
fore  Equation  (2-7)  was  used  to  compute  the  values  of  C  for  the  other 
flights.  These  values  are  listed  in  the  last  column  of  Table  2-1. 

For  each  of  the  flights  2021,  2024,  and  2042,  the  angle  0,  listed  in  the 
fourth  column  of  Table  2-1,  was  then  obtained  by  averaging  the  values  of  0 
given  by  Equations  (2-8)  and  (2-9).  (The  angle  £  computed  as  indicated  above 

0 

was  used  in  these  equations.)  The  extreme  values  of  were  obtained  by 
averaging  over  the  entire  flight  the  period  measurements  of  the  radiation 
signal  which  occurred  at  the  appropriate  times  within  the  precession  period. 

On  flights  2019  and  2027,  the  precession  and  rotation  periods  of  the 
rocket  were^  comparable  in  magnitude.  Hence,  the  method  used  above  to  deter¬ 
mine  the  variation  was  not  applicable.  The  following  method  was  used  to 
determine  the  values  of  0  for  these  flights:  Starting  at  a  time  T  ,  where 
the  counting  rate  of  the  detector  in  Channel  1  was  high,  each  value  of  the 
counts  per  pedestal  given  in  the  data  sheets  was  recorded  against  the  time  of 
occurrence  of  the  pedestal  minus  the  sum  of  Tq  and  an  integral  number  of  pre¬ 
cession  periods.  Tne  number  of  precession  periods  used  in  this  subtraction 
was  that  which  yielded  a  net  time  interval  which  was  greater  than  zero  but 
less  than  a  full  precession  period.  Also,  the  rotational  angle  was  computed 
for  ea<jh  of  these  data  points  by  multiplying  the  rotational  angular  velocity 

by  the  time  of  occurrence  of  the  pedestal,  minus  the  starting  time  T  •  A  doduc- 

o  0 

tion  of  360  from  the  rotation-angle  computations  was  made  whenever  a  complete 

rotation  was  exceeded.  After  about  300  data  points  were  recorded  in  this 

manner,  the  precession  period  was  divided  into  12  equal  intervals  of  time, 

and  the  counts  per  pedestal  which  fell  into  each  of  these  time  intervals  were 

plotted  against  the  rotation  angle  corresponding  to  that  data  point.  Each  of 
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these  graphs  therefore  yielded  the  contour  of  the  radiation  plane  (actually 
the  fold  with  the  angular  response  of  the  detector)  for  a  small  angular 
range  of  the  precession  angle.  Then  a  graph  was  prepared  giving  the  ro¬ 
tational  angle  at  which  the  peak  of  each  of  these  distributions  occurred 
versus  the  time  increment  of  the  precession  period.  This  graph  yielded 
a  sinusoidal  curve;  the  difference  between  the  maximum  and  minimum  values 
was  equal  to  26  (strictly  true  only  for  <£”  «  90°,  see  Figure  2-1). 

A  very  small  displacement  of  the  peaks  was  observed  in  these  distribution 
curves  for  the  flights  2019  and  2027,  and  the  resulting  estimates  of  the 
angles  0  are  given  in  Table  2-1:  (This  method  of  analysis  was  verified  by 
applying  it  to  flight  2024.) 

2.1.4  POTERMINATION  OF  PHASE  ANGLES  d  AND  d  . 

P  s 

The  precessional  phase  angles  d^  were  obtained  by  fitting  the  curve  given 
by  Equation  (2-6),  with  the  values  of  £  and  0  computed  as  indicated  above, 
to  the  measurements  of  the  periods  of  the  radiation  signal  within  a  precession 
period.  Such  a  curve  is  drawn  to  the  data  in  Figure  2-2. 

The  rotational  phase  angle  was  determined  for  each  group  of  data  which 
were  read-out  consecutively  (for  each  "roll"  of  data  read),  by  preparing 
graphs  of  the  counts  per  pedestal  versus  time.  These  phase  angles  were  deter¬ 
mined  in  this  manner  to  avoid  the  large  accumulative  errors. 
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2.2  ROCKET  ORIENTATION  FROM  ANTENNA  RADIATION  PATTERN 

It  was  not  necessary  to  use  the  telemetry  signal  strength  records  to  get 
the  orientation  of  the  rocket  with  respect  to  the  radiation  plane,  except  for 
the  determination  of  the  precession  period,  as  shown  by  the  above  analysis* 
These  records  were  used,  however,  to  establish  the  orientation  of  the  radia¬ 
tion  plane  in  space,  and  to  verify  the  calculations,  based  on  the  above 
methods,  of  the  angles  which  describe  the  motion  of  the  spin  axis  of  the 
vehi cle. 


Figure  2-3  depicts  a  right-handed  coordinate  system,  x,  y,  z,  whose 
origin  is  at  the  center  of  mass  of  the  vehicle.  The  x-z  plane  lies  in  a 
meridian  plane  with  the  z-axis  directed  towards  the  local  zenith  and  the 
x-axis  towards  the  south.  V  represents  a  vector  which  points  towards  a 
particular  ground  receiving  station.  Its  direction  cosines  with  respect  to 
the  x,  y,  and  z  axes  are  designated  by  the  symbols  Of  $  ^  ,  and  ^  , 

respectively.  Again  A  and  5  represent  the  constant  angular  momentum  vector 
and  the  spin  axis  of  the  vehicle;  and  £  is  the  constant  angular  displace¬ 
ment  of  3  and  A.  The  angular  displacement  of  the  S,  A  plane  from  the  node 
n(A,Z)  is  given  by  the  angle  6<Jp(t-to).  The  location  of  the  node  is  given 
by  the  angle  measured  from  the  x-axi3,  and  /(  represents  the  angular 
displacement  of  A  with  respect  to  the  z-axis. 


Now,  the  intensity  of  the  electromagnetic  radiation  recorded  by  the 
ground  station  during  the  time  interval  of  a  precession  cycle  depends  on  the 
inclination  of  the  longitudinal  axis  of  the  vehicle  with  respect  to  V.  This 
angle  will  be  represented  by  the  symbol,  x  •  14  is  given  by  the  equation, 

cosyf  =  (  0(  cos  (p  +  S  sin  (p  )  sin  cos  GO + 

(C(  sin 


(p  -  p  cos  (p)  sin sin  £(cos  ^sin  cos  n  (t"Vi 


<Kcos 
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I  Fig.  2-3  Coordinate  System  for  Determination  of  Rocket  Motion  From 

Antenna  Signal  Strength  Data 
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The  direction  cosines,  0(  ,  jS  >  ^  ,  can  be  computed  from  the  known 
geographical  locations  of  the  ground  stations  and  the  rocket.  A  relation¬ 
ship  between  5  and  K  be  obtained  as  follows:  By  correlating  the  signed 
strength  record  received  by  a  particular  ground  station  with  the  antenna 
radiation  pattern,  ^  can  be  determined  as  a  function  of  time.  Hence,  an 
integration  of  cos  "Y  (t)  over  a  precession  period  t  gives, 


+  T 


cos^f  (t)  dt  *  jjj'  Tp  cos  C  cos  /<  , 


since  the  oscillatory  terms  of  Equation  (2-10)  do  not  contribute  to  this 
integral. 

Three  additional  equations  are  required  to  evaluate  the  constants  , 
ft  ,  (jP  ,  and  tQ.  These  equations  can  easily  be  obtained  by  writing 
equations  of  the  form  of  Equation  (2-10)  for  each  of  three  ground  receiving 
stations,  utilizing  the  times  at  which  nulls  occur  (cob  =  0)  in  the 
signal  strength  records  of  these  stations. 

An  analysis  of  this  type  was  conducted  for  Test  2024.  By  utilizing 
the  magnetic  dip  and  declination  values,  at  the  surface  of  the  earth  below 
the  rocket,  given  by  the  Smithsonian  Physical  Tables,  excellent  agreement 
was  obtained  wit*'  the  results  given  by  the  analysis  described  in  Section  2.2. 
The  angles  0,  ,  and  d^  calculated  by  these  independent  methods,  were  all 

in  agreement  to  within  2  degrees.  Therefore,  at  the  positions  of  the  rockets 
in  space  which  were  used  in  these  analyses,  the  radiation  plane  was  found  to 
be  perpendicular  to  the  direction  of  the  local  geomagnetic  field. 
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3.  CALIBRATION  OF  DETECTORS 

The  geometric  functions  of  the  detectors,  G(E,  ^  yOf),  described  in 
Section  1  were  determined  by  irradiating  the  detectors  with  an  electron 
beam  of  known  energy  and  intensity  in  a  large  vacuum  chamber.  A  l80  degree 
beta- ray  spectrometer  was  used  to  obtain  a  beam  of  monoenergetic  electrons. 

With  a  SrY  source  at  the  spectrometer  entrance,  the  output  beam  energy  could 
be  adjusted  to  any  value  up  to  about  2  Mev.  The  entire  Jason  instrumentation 
package  was  placed  in  the  beam  to  take  into  account  the  scattering  of  electrons 
from  exposed  surfaces  into  detection  channels.  Although  the  electron  beam  waa 
highly  collimated,  a  combination  of  rotational  and  translational  motions  of 
each  detector  under  calibration  in  the  beam  simulated  the  effect  of  a  broad 
planar  flux  covering  the  extent  of  the  detector.  A  special  glmbal  mechanism 
was  constructed  to  support  the  instrumentation  package  at  all  possible 
spatial  orientations.  This  mechanism  facilitated  the  calibration  of  the 
detectors  as  functions  of  the  angles  'fj  and  ^  which  give  the  orientations 
of  the  detectors  with  respect  to  the  radiation  plane. 

The  radiation  plane  was  simulated  by  rotating  the  gimbal  support  mech¬ 
anism  on  a  turntable  at  constant  angular  velocity.  To  avoid  twisting  cables 
connected  to  the  geiger  tubes,  the  turntable  motor  was  stopped  and  reversed 
in  direction  each  time  the  geiger  tube  under  calibration  was  rotated  to  an 
angular  position  where  the  collimation  of  the  tube  reduced  its  counting  rate 
to  zero.  The  "thickness"  of  the  radiation  plane  was  effectively  extended  to 
cover  each  long  detector  (Anton  106C  geiger  tube)  by  moving  the  tube  through 
the  beam  along  the  long  dimension  of  the  collimator.  This  translational 
motion  was  provided  during  the  calibration  by  a  motor-driven  carriage  mounted 
on  the  gimbal  support. 

A  translation  of  the  Anton  222R  geiger-tube  detectors  through  the  beam 
was  not  required  because  the  collimator  apertures  of  these  detectors  were 
smaller  than  the  cross  sectional  area  of  the  beam. 
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Figure  3-1  illustrates  the  manner  in  which  the  orientation  angles  ^ 
and  y*'  of  the  long  detectors  were  adjusted  in  the  laboratory.  The  angle 
yj  ,  which  gives  the  displacement  of  the  spin  axis  of  the  rocket  from 
the  local  geomagnetic  line  of  force,  was  adjusted  by  inclining  the  geiger 
tube  axis  to  an  angle  7T / 2  -  ^  with  respect  to  the  horizontal  plane 
(plane  of  the  turntable).  The  angle  ,  which  denotes  the  displacement 

of  the  "central-aperture”  vector  of  the  detector  from  the  radiation  plane, 
measured  in  a  plane  perpendicular  to  the  spin  axis  of  the  rocket,  was  set 
in  the  laboratory  by  rotating  the  detector  about  the  line  L  L.  This  line 
is  shown  in  Figure  3-1.  It  is  located  in  the  entrance  plane  of  the  colli¬ 
mator,  and  it  bisects  the  aperture  along  the  long  dimension  of  the  collima¬ 
tor. 


Figure  also  depicts  the  proper  relative  locations  of  the  beam,  the 
turntable  axis,  and  the  geiger  tube.  As  shown  in  the  figure,  they  were 
arranged  such  that  the  line  L  L  of  the  detector,  the  axis  of  the  turntable, 
and  the  center-line  of  the  beam  intersected  at  the  single  point  marked  0. 
This  point  will  be  referred  to  as  the  central  irradiation  point. 

Preparation  a  calibration  run  consisted  of  adjusting  the  angles 
and  /\^>  ,  setting  the  beam  energy  E,  and  running  the  translation  carriage 
to  one  of  its  extreme  positions.  (The  detector  was  out  of  the  beam  at  the 
extreme  positions  of  the  carriage.)  The  run  was  begun  by  turning  on  a 
scaler  which  recorded  the  geiger  tube  pulses  and  a  special  electronic  unit 
which  controlled  the  rotational  and  translational  motions  of  the  detector. 
The  motions  were  fully  automatic,  with  the  carriage  advancing  a  small 
distance  S  along  line  L  L  for  each  rotation  of  the  turntable.  After 
the  carriage  traversed  the  beam  and  reached  its  other  extreme  position,  the 
detector  motions  were  stopped  and  the  scaler  was  turned  off  automatically. 
The  total  number  of  counts,  C^,  accumulated  during  the  run  was  recorded. 

To  obtain  the  background  counts,  C^,  the  run  was  repeated  with  the  beam  off 
(zero  magnet  current). 
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Schematic  Diagram  Showing  Orientation  of  Det actor  in  Electron  Beam 
The  dashed  line  LL  Is  the  axis  about  which  the  detector  le^Sefto 

Se  dXVr  m  Lrf  reSpre'  Thi»  «“  Wet.  tte  oSSSJor  of 
tne  detector  in  the  entrance  plane  as  shown 
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The  geometric  function  was  then  computed  by  means  of  the  equation, 

o(^,y,E)  a  ,  (3-1) 

where 

AS  »  linear  displacement  of  detector  per  rotation  of  turntable 
w  m  width  of  beam 

=  total  counts  during  run 

C_  *  background  counts 

D 

•*2  *1 

l(E)  =  electron  beam  intensity  (number  .  cm  •  sec  ) 

T  »  period  of  rotation  of  turntable 

This  equation  can  easily  be  derived  formally,  but  it  is  somewhat  more 
direct  and  instructive  to  verify  the  equation  purely  on  the  basis  of  physical 
considerations.  In  Section  1  the  geometric  function  of  a  detector  was 
defined  by  the  equation, 

Ri'rjfy)  =  f(E>  dE>  (3-2) 

where  R (^  ,^  )  is  the  counting  rate  of  the  detector  oriented  at  angles  Y) 

and  W  to  a  planar  flux  f(E).  The  function  G  has  units  of  cm2  when  f(E) 

'  -2-1-1 
is  given  in  units  of  number  •  cm  •  sec  .  •  Mev  and  R  is  expressed  in 

counts  •  sec  \  If  the  planar  flux  consists  essentially  of  electrons  of 

energy  Eq  within  a  small  energy  spread  AE,  then  an  integration  over  ^E 

yields 

-  G(^,y,E^)  F(E^), 

where  F(E^)  represents  the  intensity  of  the  planar  distribution  of  electrons 

O  1 

of  enevgy  Eq  within  AE  in  units  of  number  .  cm”  •  sec"  .  Moreover,  if 
the  counter  is  maintained  in  this  position  for  a  time  T,  its  counts,  C, 
would  be  given  by  the  equation 

c  =  rt  =  o(rj  eo)  F(Eq)  T.  (3-3) 
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Now,  during  the  calibration  process,  the  linear  motion  of  the  detector 

through  the  beam  effectively  widens  the  beam  to  the  extent  of  the  tube  and 

beyond.  The  total  number  of  electrons  per  unit  area  of  this  "effective" 

beam  to  which  the  detector  is  exposed  is  given  by  the  number  of  electrons 

which  traverse  a  differential  area  moving  through  the  width,  w,  of  the 

actual  beam  at  the  rate  of  (  4  6  cos  Yj  )/T.  (See  Figure  3-1)*  Since  the 

beam  intensity  of  the  actual  beam  is  I(Eq),  this  number  is  I  •  w  •  (T/^S» 

cos  “Yj  ).  But  since  the  tube  was  rotated  during  the  calibration,  this  areal 

density  of  electrons  was  uniformly  distributed  in  angle  around  the  detector. 

It  is  therefore  equivalent  to  an  irradiation  of  the  detector  with  a  planar 

-2  -1 

flux  of  F  electrons  •  cm  .  sec  for  a  time  duration  T.  Or, 

FT  =  I  w  T/(  A  S  cos  'Yj  ).  (3-4) 

When  this  value  of  F  is  put  into  Equation  {3-3) >  the  validity  of  Equation 
(3-1)  is  verified. 

Since  the  apertures  of  an  Anton  222R  detectors  were  much  smaller  than 
the  cross  sectional  area  of  the  beam,  each  of  these  detectors  was  calibrated 
by  locating  the  center  of  the  entrance  aperture  of  the  collimator  at  the 
central  irradiation  point  and  rotating  the  detector  turntable  while  scaling 
the  counts.  The  geometric  function  was  then  computed  by  using  Equation  (3-3) 
in  the  form, 

0('7’Y’E)  -  -£(E)tB  -  (3-5) 

where  N  is  the  number  of  times  the  detector  was  rotated  in  the  beam  to  obtain 
good  statistics  (N  100). 
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3.1  EXPERIMENTAL  APPARATUS 

The  experimental  equipment  used  in  the  calibration  of  the  detectors, 
except  for  a  rack  of  electronics  equipment,  is  shown  diagrammatically  in 
Figure  3-2  and  pictorially  in  Figure  3-3.  The  equipment  can  be  grouped 
into  three  principal  categories:  the  beta- ray  spectrometer;  the  planar-beam 
simulation  platform,  which  5ncludes  the  gimbal  support  mechanism  as  well  as 
the  turntable  and  translation  carriage  mechanisms;  and  the  electronics  equip¬ 
ment. 


3.1.1  THE  BETA-RAY  SPECTROMETER 

The  beta-ray  spectrometer  is  shown  at  the  left  of  the  chamber.  This 
consisted  of  a  180-degree  magnetic  deflection  spectrometer  using  a  one  milli- 
curie  Sr-Y  source.  The  spectrometer  was  mounted  on  an  adjustable  platform 
and  located  on  the  small  table  shown  in  Figure  3-2.  The  beta-decay  source 
and  the  collimation  slits  which  selected  the  radii  of  curvature  (  p  pc.  3  cm) 
of  the  electrons  were  mounted  in  an  aluminum  vacuum  tight  housing  of  outside 
dimensions  2-1/2"  x  3-1/2"  x  l/2",  and  this  "wafer"  was  inserted  into  the 
gap  of  the  magnet.  A  silphon  link  was  used  to  connect  the  wafer  to  the 
vacuum  tank.  In  order  to  rapidly  suppress  the  stray  field  of  the  magnet 
along  the  path  of  the  electron  beam,  a  one-inch  thick  cylinder  of  cold-roll 
iron  was  sealed  between  the  wafer  housing  and  the  silphon,  as  shown  in 
Figure  3-2.  A  window  was  cut  through  the  iron  shield  to  enable  the  electron 
beam  to  enter  the  vacuum  chamber.  The  flexible  link  and  the  adjustable 
magnet-support  platform  facilitated  the  positioning  of  the  beam  in  the  tank. 

After  the  beam  entered  the  vacuum  chamber,  it  was  further  collimated  to 
limit  its  cross  section  to  3"  x  9/l6"  at  the  central  irradiation  point  "0" 
shown  in  Figure  3-1.  A  schematic  diagram  of  the  spectrometer  and  the  beam 
is  shown  in  Figure  3-4.  The  spatial  and  energy  distributions  of  the  beam 
are  described  in  Section  3.2 

A  special  variable-current  power  supply  was  constructed  for  the  magnet 
of  the  spectrometer.  This  supply  delivered  up  to  7  amps  at  2k  volts,  which 
weis  more  than  sufficient  to  reach  the  end-point  energy  of  the  SrY  source. 

The  desired  magnet  current  was  adjusted  by  reading  the  potential  across  an 
80-mil  ni chrome-wire  shunt  with  a  Rubicon  portable  potentiometer. 
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Fig.  3-2  Diagram  of  Experimental  Setup 
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3.1.2  PLANAR-BEAM  SIMULATION  PLATFORM 

In  this  section  a  description  is  Given  of  the  "target"  platform  which 
is  depicted  at  the  center  of  the  vacuum  chamber  in  Figure  3-2.  This  plat¬ 
form  serves  a  dual  function:  it  supports  the  detectors  at  the  required 
orientations  with  respect  to  the  simulated  radiation  plane,  and  it  rotates 
and  translates  the  detectors  in  the  electron  beam  to  effect  an  irradiation 
of  the  detectors  by  a  planar  flux.  A  diagram  of  this  assembly  is  shown  in 
Figure  3-5. 

The  turntable  was  located  at  the  bottom  of  the  cylindrical  "well"  of 
the  vacuum  chamber  as  indicated  in  Figure  3-2.  An  O-ring  seal  was  used  to 
feed  the  shaft  of  the  turntable  through  the  base  plate  of  the  chamber. 

The  turntable  drive  mechanism  together  with  a  cam-microswitch  arrangement 
which  controlled  the  maximum  angle  of  rotation  of  the  turntable  in  either 
direction  were  located  beneath  the  chamber.  A  Bodine  synchronous-capacitor 
motor  drove  the  turntable  through  reducing  gears  at  the  uniform  rate  of 
0.25  rpm.  This  equipment  under  the  chamber  is  shown  in  Figure  3-6. 

The  gimbal  mechanism,  which  rides  on  the  turntable  and  supports  the 
translation  carriage  to  which  the  instrumentation  package  is  attached,  is 
illustrated  in  Figure  3-5*  Also,  pictures  of  this  device  are  shown  in 
Figures  3-7a  and  3-7h.  Figure  3-7a  affords  a  good  view  of  the  motor  and 
screw-drive  mechanism  which  move  the  detector  package  across  the  electron 
beam.  A  study  of  Figures  3-5>  3-7a>  and  3-7^  reveals  the  manner  in  which 
the  detector-orientation  adjustments  were  made.  The  angle  7^  was  set  by 
adjusting  the  positions  of  the  clamps  on  the  two  C-bars  mounted  to  the  base 
plate.  The  angle  was  adjusted  by  rotating  the  carriage  mount  structure 

about  a  shaft  fixed  to  the  brace  which  rides  on  the  C-bars. 

The  carriage -drive  motor  was  a  Bodine  two-phase  servomotor  which  hnd  on 
output  of  22  revolutions/minutc  through  a  gear-reducer  system.  Since  it  un¬ 
necessary  to  operate  this  motor  in  the  vacuum  system,  special  precaution:, 
were  taken:  the  original  oil  in  the  motor  was  replaced  with  Silicon  gr 
and  a  heat  sink  for  the  motor  was  constructed.  The  heat  sink  provision  con¬ 
sisted  of  constructing  good  thermal  paths  for  the  conduction  of  the  motor  heat 
to  the  structure  of  the  gimbal  mechanism. 
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Fig.  3-6  Turntable  Drive  Mechanism  and  Limit  Microswitches  Under  Vacuum  Chamber 
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Fig.  3- 7a  Gi ruble  Support  Showing  Carriage-Drive  Mechanism 
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Fig.  3-7b  Additional  View  of  Gimble  Support 
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3.1.3  ELECTRONICS  EQUIPMENT 

The  automatic  operation  of  the  rotation-translation  motions  of  the 
detectors  was  controlled  by  a  master  control  unit.  The  operation  was  such 
that  the  carriage  was  moved  in  preselected  movements  of  l/l6-inch,  l/Q-inch, 
or  l/4-inch  every  time  a  turntable  limit  switch  wan  hit  by  the  pre-set  cams 
on  the  turntable  axis.  A  block  diagram  of  the  logic  of  the  automation 
system  which  performed  this  function  is  shown  in  Figure  3 -$»  and  the 
operation  of  the  equipment  during  a  calibration  run  is  outlined  below. 

(1)  The  cams  on  the  axis  of  the  turntable  were  a/1  justed  for 
the  desired  limit  of  rotation  of  the  turntable  in  each 
direction.  These  limits  were  set  by  observing  the  rota¬ 
tion  angles  at  which  the  counting  rate  of  the  detector 
under  calibration  dropped  to  its  cosmic-ray  background 
level. 

(2)  The  "Electrical  Override"  on  the  master  control  unit 
was  then  operated  to  move  the  carriage  to  the  starting 
position.  At  this  position  the  detector  was  just  out  of 
the  beam. 

(3)  The  translation  movement  of  the  carriage  (l/l6",  1/8",  or 
1/4")  for  the  run  was  selected  by  adjusting  the  "Stepping 
Device  Preselector"  switch. 

(4)  The  run  was  then  started.  The  turntable  rotated  in  a 
certain  direction  until  a  limit  switch  was  hit. 

(5)  The  actuation  of  the  limit  switch  initiated  the  following 
sequence  of  events: 

•  The  turntable  motor  was  stopped. 

•  The  stepping  switch  was  energized,  and  it 

applied  a  step  voltage  to  the  input  of  a 
servo-amplifier.  The  magnitude  of  this 
voltage  controlled  the  incremental  motion 
of  the  carriage. 
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•  The  output  of  the  servo-amplifier  energized 
the  carriage  motor  until  the  new  position 
Of  the  carriage  was  reached. 

•  The  correct  position  of  the  carriage  was 
sensed  by  the  balancing  of  the  voltage 
applied  to  the  servo-amplifier  by  the 
stepping  switch  and  a  voltage  controlled  by 
a  micropot  connected  to  the  screw-drive 
shaft  of  the  carriage. 

•  After  the  correct  carriage  position  was 
reached,  the  turntable  motor  was  again 
actuated,  but  with  reversed  polarity, 
effecting  the  rotation  of  the  detector  in 
the  opposite  direction. 

(6)  This  process  was  repeated  each  time  a  limit  switch  was 
actuated  until  the  carriage  reached  its  opposite  limit 
position  across  the  beam. 

(7)  After  this  position  was  reached,  the  detector  motions  were 
stopped,  and  the  scaler  which  recorded  the  total  number  of 
detector  pulses  during  the  run  was  disengaged. 

Thirty  to  seventy  rotations  through  the  beam  were  made  during  a  run.  The  time 
duration  of  a  run  varied  from  30  to  60  minutes. 
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3.2  ALIGIWENT  AND  CALIBRATION  OF  ELECTRON  BEAM 

In  preparing  the  spectrometer  for  the  detector  calibration  work,  the 
fan  plane  of  the  beam  (the  medium  plane  through  the  spectrometer  magnet 
gap)  was  first  aligned  with  the  center-line  of  the  turntable  axis.  Tills 
alignment  was  accomplished  by  adjusting  the  position  of  the  spectrometer 
outside  the  chamber  while  monitoring  the  location  and  distribution  of  the 
beam  inside  the  chamber. 

Energy  spectra  of  the  electrons  were  measured  throughout  the  cross 
section  of  the  beam  by  means  of  a  plastJ  c  scintillation  spectrometer. 

Pulse  height  distributions  from  the  spectrometer  were  recorded  by  a  100- 
channel  pulse  height  analyzer.  Spatial  distributions  of  th^  intensity  as 
well  as  the  spectra  of  the  electrons  in  the  beam  cross  section  were  measured 
by  mounting  a  collimating  slit  in  front  of  the  plastic  scintillator  and 
moving  the  spectrometer  through  the  beam  with  the  translation  carriage 
mechanism.  By  setting  the  gimbal  settings  to  Tj  =  90°  and  Aj/'  =  0, 
horizontal  cuts  through  the  bean  were  taken.  Such  traversals  were  made  at 
various  elevations  with  respect  to  the  central  irradiation  point  to  test 
the  vertical  alignment  of  the  fan  plane  of  the  beam.  With  the  gimbal  settings 
at  ^  =  0  and  /\j/s  =  0,  vertical  cuts  through  the  beam  were  taken. 

With  no  external  collimation  of  the  beam,  the  results  of  these  measure¬ 
ments  revealed  the  following  information: 

(1)  The  edges  of  the  beam  were  quite  sharp,  beinp  "smeared" 
over  a  distance  which  was  well  within  the  width  of  the 
detector  slit  (3/l6"). 

(2)  The  cross  section  of  the  beam  had  a  width  of  about  5/0"  and 
a  height  of  6". 

(3)  The  enerGy  spectrum  was  quite  uniform  over  the  entire  cross 
section:  the  Greatest  departure  from  uniformity  was 
observed  at  the  lowest  electron  energy  of  interest  (E 

150  kev);  here,  the  displacement  of  the  peaks  of  the  pulse- 
heiGht  spectra  over  the  6"  heip,ht  of  the  beam  was  about  5$, 
and 
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(4)  The  intensity  distribution  along  the  height  of  the  beam 
was  far  from  being  uniform:  it  varied  rather  linearly 
along  the  height,  differing  by  nearly  a  factor  of  two 
from  one  end  to  the  other. 


The  non-uniform  vertical  distribution  of  the  electron  beam  intensity 
was  unfortunate,  and  it  was  necessary  to  restrict  the  height  of  the  beam 
to  a  region  where  the  linear  approximation  of  the  intensity  variation  was 
good.  A  linear  variation  was  not  considered  to  be  troublesome  since  the 
detectors  were  symmetrical  above  and  below  the  central  irradiation  point 
"O'. 

Since  it  was  necessary  for  the  beam  to  cover  the  entire  aperture  of 
the  detector  in  the  vertical  direction,  the  minimum  height  of  the  beam  was 
dependent  on  the  orientation  angle  and  on  the  width  of  the  beam.  This 
height  was  found  to  be  about  1-1/4"  for  the  maximum  range  in  “Yj  (69°  to 
lll^ given  by  the  analysis  of  the  body  motions  of  the  rockets  and  for  a 
beam  width  of  9/16".  Hence,  the  external  collimator  shown  in  Figure  3-4 
was  arranged  to  restrict  the  width  of  the  beam  to  9/l6"  and  the  height  of 
the  beam  to  about  +1"  from  the  central  position. 

The  horizontal  distribution  of  the  beam  across  the  central  position 
was  then  measured,  and  the  result  is  shown  in  Figure  3-9*  The  solid-line 
trapezoid  drawn  in  the  figure  depicts  the  fold  of  a  3/16"  slit  with  a  9/16- 
wide  pedestal-shaped  beam  distribution.  The  vertical  distribution  of  the 
beam  is  shown  in  Figure  3-10.  The  origin,  point  0,  of  the  abscissa  corres¬ 
ponds  to  the  central  irradiation  point. 

After  the  beam  was  aligned  and  calibrated,  the  intensity  and  energy  of 
the  beam  versus  magnet  current  were  measured  with  a  plastic  scintillation 
spectrometer.  A  collimator  with  an  aperture  3/4"  wide  and  1-1/2  high  was 
placed  in  front  of  the  2"  diameter  plastic  scintillator,  and  the  spectrometer 
was  mounted  on  the  gimbal  mechanism  such  that  the  center  of  the  aperture 
coincided  with  the  central  irradiation  point  "0"  shown  in  Figure  3-1  and 
Figure  3-5.  In  this  condition,  therefore,  the  entire  width  of  the  beam  and 
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Fig.  3-9  Horizontal  Intensity  Distribution  of  Beam  ■ 
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the  height  of  the  beam  from  -3/V'  to  +3/V'  shown  In  Figure  3-10  fell  within 
the  aperture  of  the  spectrometer.  Pulse-height  spectra  were  then  obtained 
with  a  100-channel  analyzer  as  a  function  of  magnet  current.  For  each 
value  of  the  magnet  current  used,  the  beam  energy  was  taken  to  be  that 
given  by  the  peak  of  the  pulse-height  spectrum  after  the  background  spectrum 
was  subtracted.  Various  ^-emitting  radioactive  sources  were  used  to 
calibrate  the  scintillation  spectrometer.  The  intensity  of  the  beam  was 
computed  from  the  total  number  of  counts  in  the  distribution  above  the  back¬ 
ground.  The  energy  of  the  beam  as  a  function  of  the  potentiometer  reading 
of  the  voltage  across  a  standard  shunt  is  given  in  Figure  3-H*  The  inten¬ 
sity  of  the  beam,  also  as  a  function  of  the  potentiometer  reading,  is  given 
in  Figure  3-12. 

By  cycling  the  magnet  current  before  adjusting  the  current  to  the 
desired  value,  the  beam  energy  was  found  to  be  reproducible  to  about  2$. 

A  recalibration  of  the  beam  after  a  month  of  operation  revealed  a  negligible 
change  in  the  beam  characteristics. 

The  beam  energy  as  a  function  of  magnet  current  was  also  computed  by 
measuring  the  magnetic  field  in  the  center  of  the  magnet  gap  with  a  Vjo 
flip-coil  magnetometer  and  the  mean  radius  of  curvature  allowed  by  the 
spectrometer  stops.  These  computations  agreed  well  with  the  spectrometer 
measurements  when  a  radius  of  2.90  cm  was  used  rather  than  the  measured 
value  of  3*00  cm.  This  error  probably  resulted  from  the  fact  that  the  peak 
field  intensity  rather  than  the  mean  intensity  along  the  path  of  the 
electrons  was  used.  Further  agreement  with  the  energy  calibration  of  the 
beam  was  observed  by  preparing  a  Fermi  plot  of  the  intensity  distribution 
data.  This  plot,  shown  in  Figure  3-13,  is  seen  to  define  a  straight  line 
which  yields  the  end-point  energy  of  the  SrY  source  used  in  the  spectrometer 
to  within  a  few  percent. 

3.3  DETECTOR  CALIBRATION  RESULTS 

The  method  used  in  the  analysis  of  the  data  obtained  with  the  Wallops- 
Island  launched  rockets  required  knowledge  of  the  geometric  functions  of 
the  detectors  in  the  orientation  range  69°^  ^  111°  and  -20°^  20C. 
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Fig.  3-11  Electron  Beam  Energy  as  Function  of  Magnet  Current.  Current 
reading  is  in  terms  of  the  potentiometer  reading 
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Fig.  3-12  Electro  Beam  Intensity  Over  1.  5-in.  Beam  Height 
(Elec.rons/Sec)  as  a  Function  of  Magnet  Current 
(Potentiometer  Reading) 
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Fig.  3-13  Fermi  Plot  Given  by  Beta-Ray  Spectrometer  Data.  Abscissa  gives  total 
energy  In  units  of  rest  mass  energy  of  the  electron 
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The  range  in  O^j  was  given  by  the  determination  of  the  body  motions  of  the 
rockets  in  free  flight.  The  range  in  of  +  20°  was  important  because 

only  those  counting  rates  of  the  detectors  which  were  read  out  in  this 
range  of  Oj/'  were  used  in  the  analysis.  Counting  rates  which  were 
sampled  in  flight  while  the  detectors  were  oriented  at  larger  angles  of 
were  not  used  becaise  it  was  believed  that  large  errors  would  have  resulted 
from  the  non-planar  distribution  of  the  trapped  electrons. 

The  geometric  functions  measured  in  the  laboratory  were  obtained  for 
various  values  of  the  parameters  'Vj  ,  'tj/'  ,  and  E.  The  data  below  reveal, 
however,  that  the  geometric  functions  of  the  long  tubes  were  quite  insensi¬ 
tive  (less  than  5$)  to  the  orientation  angle  7J  within  the  range  given 
above.  This  dependence  was  therefore  neglected.  Moreover,  since  the  count¬ 
ing-rate  samples  which  occurred  within  the  0^  interval  of  +20°  were  averag'  d 
over  a  number  of  spin  cycles  of  the  rocket,  the  geometric  functions  used  to 
unfold  the  data  were  also  averaged  over  this  interval  of  in  a  manner 
discussed  below.  The  end  result,  therefore,  was  a  geometric  function  for 
each  detector  which  was  dependent  on  the  electron  energy  alone. 


The  counting  rates  of  the  small  detectors  were  analyzed  in  a  somewhat 
different  manner,  as  discussed  below.  But  that  method  also  required  an 
average  of  the  geometric  functions  over  the  interval  -20°^  ^  20°. 

The  calibration  data  are  given  below,  with  the  results  for  the  long 
detectors  preceding  those  for  the  small  detectors.  Within  each  of  those 
groups  of  detectors,  the  detectors  are  considered  in  the  order  of  their 
increasing  threshold  energies. 


3.3.1  CALIBRATION  OF  CHANNELS  1  AND  5 

Hie  detectors  in  Channels  1  and  5  were  identical  in  construction.  Each 

O 

consisted  of  an  Anton  106C  Geiger  tube  (30  mg/cm  -  steel  wall  thickness) 
with  a  brass  collimator  forming  six  rectangular  windows  along  the  length  of 
the  tube  (see  diagram  in  Figure  1-1 ). 

First  the  orientation  angle  of  the  detector  in  Channel  1  was  set 

equal  to  90°  (longitudinal  axis  of  detector  in  radiation  plane),  and  calibra¬ 
tion  runs  were  made  at  a  number  of  values  of  beam  energy  for  each  of  the 
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following  settings  of  :  -25.2°,  -19.8°,  -12.6°,  -3«6^,  +3*8°,  +10.8° 

and  +19.8°.  Using  Equation  (3-1)  the  geometric  functions  were  computed 
from  the  data  obtained  at  each  orientation  of  the  detector,  and  the 
functions  were  plotted  against  electron  energy.  A  graph  of  the  results 
for  the  position  'Yj  -  90°,  =  3.6°  is  shown  in  Figure  3-14.  The  error 

bars  drawn  at  the  data  points  are  standard  deviations  due  to  counting 
statistics.  Error  bars  are  not  shown  at  the  lower  energy  points  where  the 
errors  are  smaller  than  the  circles  depicting  the  data.  The  threshold 
energy  of  this  detector  is  seen  to  lie  at  about  200  kev. 


To  obtain  an  average  of  these  geometric  functions  over  the  angle  /y/'  , 
the ' functions  were  read  off  the  curves  drawn  through  the  data  as  shown  in 
Figure  3-lU  at  various  values  of  electron  energy  and  plotted  against  the 
angle  .  A  graph  of  the  resulting  curves,  for  a  few  of  the  energy  para¬ 
meters  used,  is  shown  in  Figure  3-15*  The  geometric  function  for  each  of 
the  energy  parameters  used,  was  then  averaged  over  the  interval  -20°  £ 

20°  by  means  of  a  planimeter.  The  geometric  function  given  by  this  averaging 
process  is  shown  in  Figure  3-l8. 

The  result  of  a  cali oration  run  taken  at  OYj  =  72°,  =  0  is  shown  in 

Figure  3 -17.  Within  the  statistics  of  the  measurements,  the  geometric 
function  given  in  this  figure  is  indistinguishable  from  that  given  by  extra¬ 
polating  the  measurements  taken  at  'Yj  =  90°,  ^  =-3.6°  (see  Figure  3-1*0 

to  the  peak  response  at  /\jf'  =  0  (see  Figure  3-15). 

Two  additional  calibration  runs  were  taken,  one  at  Yj  =  72°,  /|^  =  0 
and  the  other  at  Yj  =  72°,  =  19.8°.  During  these  runs,  three  of  the 

six  collimator  windows  of  the  detector  (see  Figure  1-1)  were  covered  with 
l/8"  of  copper,  as  was  done  on  some  of  the  flights.  The  covered  windows 
were  the  window^  at  the  extreme  end  of  the  detector  (end  towards  nose  of 
rocket)  and  the  other  two  were  the  adjacent  windows  at  the  opposite  end  of 
the  detector.  The  ratio  of  the  geometruc  functions  of  the  plugged  and 
unplugged  detectors  measured  at  the  same  orientation  to  the  beam  (  Yj  =  72°, 
a  0),  was  found  to  be  0.652.  Moreover,  this  value  agreed  with  the  ratio 
of  the  geometric  function  of  the  plugged  detectors  at  Yj  a  72°,  ^  *  19.8° 

to  that  of  the  unplugged  detector  at  Yj  =  90°,  W'  a  19.8°,  yielding 
further  evidence  of  the  insensitivity  of  the  geometric  functions  to  relati¬ 
vely  small  departures  of  Yj  from  90°. 
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Table  1  in  the  Appendix  lists  the  energy-dependent  geometric  function 
of  the  plugged  detector  in  Channel  1.  It  is  the  geometric  function  shown 
in  Figure  3“l6  multiplied  by  O.652. 

3.3.2  CALIBRATION  OF  CHANNEL  7 

The  detector  in  Channel  7  also  consisted  of  an  Anton  106C  geiger  tube, 

but  this  tube  was  not  deliberately  collimated.  It  was  shielded  with  a 

2 

cylindrical  sheath  of  aluminum  about  150  mg/can  thick. 

Visual  observations  revealed  that,  throughout  the  orientation  range 
-69°  £  7j  ^  111°,  -20  £■  ty1  £  20°  and  beyond,  the  aperture  of  the 
detector  in  a  planar  flux  was  unchanged.  Measurements  of  the  counting  rate 
of  the  tube  at  various  electron  energies  revealed  that  this  was  indeed  the 
case.  For  *7/  =  90°  these  counting  rates  were  found  to  be  constant  over  a 
total  range  of  ^  of  about  60°.  Hence,  an  averaging  of  the  geometric  function 
of  this  detector  over  was  not  necessary. 

The  energy-dependent  geometric  function  of  this  detector  is  shown  in 
Figure  3-l8a.  The  threshold  energy  is  about  470  kev.  Since  these  data  In¬ 
dicate  that  the  response  of  the  tube  is  still  rising  appreciably  at  the 
highest  electron  beam  energy  used,  the  cylindrical  shield  was  removed  to 
measure  the  geometric  function  of  the  "bare"  tube.  This  bare  tube  value  was 
considered  to  be  a  reasonable  value  which  the  shielded  detector  would  approach 
at  higher  energies.  Hence,  the  curve  drawn  to  the  data  was  extrapolated  to 
the  "bare"  tube  value  as  shown  in  Figure  3-l8b.  This  geometric  factor  is  also 
listed  against  energy  in  Table  A-l. 

3.3.3  CALIBRATION  OF  CHANNEL  3 

The  detector  in  Channel  3  again  was  an  Anton  106C  geiger  tube,  not 
deliberately  collimated,  but  encased  within  a  cylindrical  shield  of  aluminum, 

o 

kOO  mg/ cm  in  thickness. 

Again,  visual  as  well  as  experimental  observations  revealed  that  the  geo¬ 
metric  function  of  this  detector  should  be  constant  in  a  planar  flux  over 
the  range  of  orientation  of  the  detectors  given  above.  The  experimental 
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results  obtained  at  two  different  values  of  Yj  but  with  ^  anywhere  in 
the  interval  -25°^  25°  are  shown  in  Figure  3-19*  The  results  indi¬ 

cated  by  the  circles  were  obtained  at  Yj  ■  90°,  and  the  crossed  lines 
show  the  results  of  the  measurements  at  Yj  a  720.  The  threshold  energy 
of  this  detector  is  at  9^0  kev. 

Unfortunately,  the  limited  beam  energy  of  the  spectrometer  highly 
restricts  the  full  calibration  of  this  detector.  It  was  necessary  to  extra¬ 
polate  the  geometric  function  of  this  detector  to  that  of  the  "bare"  tube 
value.  Since  the  extrapolation  was  extreme,  a  study  was  undertaken  to  Bcek 
a  physical  basis  for  the  extrapolation. 

The  following  empirical  method  was  devised  for  this  purpose.  (It  was 
also  used  to  compute  the  geometric  function  of  the  detector  in  Channel  6.) 

The  measured  values  of  the  geometric  functions  of  the  detectors  in  Channels 
7  and  3  were  divided  by  the  "bare"  tube  value  to  obtain  the  effective  trans¬ 
mission  function  of  each  detector.  These  transmission  curves  were  found  to 
be  systematically  lower  than  the  transmission  functions  for  normal -incidence 
electrons  on  plane  absorbers  of  corresponding  materials  and  thicknesses. 

(The  normal-incidence  transmission  data  of  Marshall  and  Ward^  and  Hereford 
L  . 

and  Swann  were  used  in  this  analysis.)  The  diminution  of  the  tube  trans¬ 
mission  from  the  normal  transmission  values  at  each  electron  energy  was 
caused  by  the  particular  distribution  of  the  incidence  angles  of  the  electrons 
on  the  tube  resulting  from  the  irradiation  of  the  cylindrical  detector  by  a 
planar  flux.  Hence,  the  ratio  of  the  tube  transmission  to  the  normal  trans¬ 
mission  plotted  against  energy  should  give  a  cufvc  which  reflects  the  combinod 
effects  of  the  shape  of  the  detector  and  the  angular  distribution  of  the 
flux.  If  this  is  true,  then  detectors  with  different  shield  thicknesses  but 
of  the  same  shape  and  exposed  to  the  some  flux  should  yield  identical  ratios 
of  tube  transmission  to  noimal-incidence  transmission  as  functions  of 
electron  range  in  units  of  absorber  thickness.  Figure  3-20  shows  that  the 
agreement  of  these  ratios  obtained  from  the  calibration  data  of  the  detectors 
in  Channels  7  and  3  is  fairly  good.  The  solid  curve  drawn  through  the  data 


3.  Marshall  and  Ward;  Canadian  Journal  of  Research;  A  1 30  (1937) 

4.  Hereford  and  Swann;  Fhys.  Rev.  'J6,  727  (1950 ) 
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given  by  both  counters  was  therefore  used  to  extrapolate  the  geometric 
function  of  the  detector  in  Channel  3  to  higher  electron  energies. 

The  final  energy-dependent  geometric  function  of  this  detector  is 
given  in  Table  Awl. 

3.3.4  CALIBRATION  OF  CHANNEL  6 

The  detector  in  Channel  6  also  consisted  of  a  shielded  Anton  I06C  geiger 
tube  with  no  deliberate  collimation.  This  counter  was  shielded  with  n 

p 

cylindrical  absorber  of  brass  2000  mg/ ctn  thick. 

In  this  case  the  threshold  energy  of  the  detector  (4.3  Mev)  precluded 
an  experimental  evaluation  of  any  portion  of  the  geometric  function  of  the 
detector  with  the  beta-ray  spectrometer.  Hence,  the  geometric  function  was 
computed  by  means  of  the  empirical  method  described  above.  For  this  purpose 
the  curve  of  Figure  3-20  was  slightly  modified  by  plotting  the  transmission 
ratios  against  electron  range  in  copper  in  units  of  the  copper-range  of 
electrons  at  the  threshold  energies  of  the  detectors  in  Channels  7  and  3* 

This  curve,  together  with  the  Hereford  and  Swann  normal-incidence  transmission 
data  for  the  absorber  (copper)  thickness  used  in  Channel  5,  were  used  to 
compute  the  transmission  function  of  that  detector.  The  geometric  function 
was  then  obtained  by  multiplying  this  transmission  function  by  the  "bare" 
tube  value  of  the  detector  in  Channel  7» 

The  resulting  geometric  function  is  shown  in  Figure  3-21  and  listed  in 
Table  Awl. 


3.3.5  CALIBRATION  OF  CHANNEL  2 

The  detector  in  Channel  2  consisted  of  an  end-window  Anton  222R  geiger 

tube,  with  a  flat  aluminum  absorber,  about  20  mg/cm  thick,  placed  over  the 

window  of  the  detector,  the  detector' was  collimated  by  means  of  a  circular 
2 

hole,  about  0.10  cm  in  area,  drilled  through  a  brass  "plug"  at  an  angle  of 
30°  with  respect  to  the  axis  of  the  geiger  tube  (see  Figure  1-1 ). 
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Fig.  3-21  Geometric  Function  of  Detector  in  Channel  6 
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Since  the  collimator  aperture  was  circular,  the  geometric  function  of 
this  detector  was  dependent  only  on  the  an^ailar  displacement  of  the  colli¬ 
mator  axis  from  the  radiation  plane .  The  data  fram  this  detector  were 
analyzed  by  using  the  counting-rate  samples  taken  when  this  angular  dis¬ 
placement  was  less  than  +10°  as  discussed  in  Section  *<-.  This  10°  displace¬ 
ment  of  the  axis  corresponds  closely  to  'Tjf  ~  -£0°  for  7j  =  TT/2.  Hence, 

the  geometric  function  used  for  the  analysis  of  the  data  was  that  measured 
for  'yj  «  l\/2  and  averaged  over  -20° £  20°.  During  portions  of  the 

precession  cycle  of  a  rocket  where  Yj  ^  90°,  this  method  of  averaging  the 
geometric  function  does  not  give  the  correct  average  of  the  function  over 
the  displacement  an&Le;  however,  the  maximum  deviation  due  to  this  cause 
is  6$,  small  in  comparison  with  the  uncertainty  in  the  collimation  of  the 
tube  discussed  below. 


The  detector  was  calibrated  by  setting  Yj  =  90°  (longitudinal  axis  of 
geiger  tube  in  radiation  plane)  and  rotating  the  detector  in  the  beam  with 
the  entrance  aperture  of  the  collimator  centered  at  the  central  irradiation 
point*  At  each  beam  energy,  a  calibration  run  consisted  of  about  100  rota¬ 
tions  of  the  turntable*  The  orientation  angle  was  adjusted  by  rotating 

the  detector  about  an  axis  which  was  in  the  radiation  plane  and  which  went 
through  the  central  irradiation  point*  Calibration  runs  were  taken  for  /\f^ 
ancles  of  -51.0°,  -36.5°,  -22.6°,  -10.4°,  0.0°,  8.8°,  18.2°,  32. 6°,  47.4°, 
and  6l.3°.  Equation  3-5  was  used  to  compute  the  geometric  functions.  The 
measurements  of  the  geometric  functions  at  Yj  =  90°,  =  18.2°  is  shown 

in  Figure  3-22.  The  threshold  energy  is  seen  to  fall  at  145  kev. 


The  statistical  accuracy  of  the  data,  as  shown  in  Figure  3-22,  was  con¬ 
siderably  poorer  than  that  of  the  long-detectors.  Therefore,  instead  of 
averaging  the  geometric  functions  over  the  angle  Off  as  was  done  for  Channel 
1,  namely,  by  plotting  the  functions  against  using  energy  as  a  parameter, 

all  of  the  data  points  were  used  by  plotting  the  integral  over  energy  of  the 
geometric  functions  against  .  A  graph  of  the  integrated  geometric 
functions  versus  is  shown  in  Figure  3-23*  The  average  value  ol'  this 

response  function  over  the  interval;  -20°^  20°  divided  by  the  response 

at  the  peak  is  0.932.  Hence,  by  multiplying  the  geometric  functions  obtained 
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at  ^  s  90°,  Ajf'  m  0°  by  this  ratio,  the  final  averaged  geometric  function 
for  the  detector  was  obtained.  This  is  shown  in  Figure  3*24  and  listed  in 
Table  A-l. 

Unfortunately,  the  circular  aperture  of  this  detector,  as  well  as  those 
of  the  detectors  in  Channels  8  and  4,  was  partially  obstructed  by  a  titanium 
shield  which  was  placed  over  the  geiger  tube  package.  The  magnitude  of  ^he 
obstruction  depended  on  the  orientation  of  the  counter  in  the  package. 
Accurate  information  on  the  positional  setting  of  counter  is  not  available. 
An  uncertainty  in  the  geometric  function  of  this  tube  of  15$  has  been  esti¬ 
mated  on  the  basis  of  the  probable  variation  of  the  position  of  the  counter. 


3.3.6  CALIBRATION  OF  CHANNELS  8  AND  4 

The  detectors  in  Channels  8  and  4  were  identical  in  construction  to 

that  in  Channel  2,  but  the  absorbers  placed  in  front  of  the  Anton  222R 

9  2 
windows  differed.  The  detector  in  Channel  8  had  an  aluminum  absorber  150  mg/cm 

thick,  and  the  detector  in  Channel  4  had  a  brass  absorber,  400  mg/cn  in 

thickness. 

Since  the  geometry  of  each  of  these  detectors  was  the  same  as  that  of 
the  detector  in  Channel  2,  only  the  transmission  factors  of  these  detectors 
were  expected  to  vary.  Hence,  the  transmission  functions  alone  were  measured 
for  these  detectors. 

The  transmission  function  of  each  detector  was  measured  by  aligning  the 

collimator  in  the  beam  at  the  central  irradiation  point  and  dividing  the 

resulting  counting  rate  of  the  detector,  above  cosmic  ray  background  by  the 

product  of  the  beam  intensity  and  the  area  of  the  aperture.  Figure  3-25 

shows  the  results  obtained  for  the  detector  in  Channel  8.  These  data  agree 

quite  well  with  the  Marshall  and  Ward  measurements  (transmission  of  electrons 

at  normal  incidence  through  150  mg/cm  of  aluminum),  which  are  also  given  in 

the  figure.  The  threshold  energy  is  400  kev.  The  transmission  measurements 
o 

with  the  400  mg/cm  brass  absorber  together  with  the  transmission  curve  given 
by  the  Marshall  and  Ward  data  for  normal-incidence  electrons  through  400  mg/cm 
of  aluminum  are  shown  in  Figure  3-26.  These  data  agree  as  well  as  can  be 
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expected  in  view  of  the  different  absorber  materials  used.  Since  a  rela¬ 
tively  short  extrapolation  of  the  data  to  higher  energies  was  required,  the 
Marshall  and  Ward  curve  was  used  to  compute  the  geometric  functions  for 
energies  greater  than  2.0  Mev.  The  threshold  energy  of  this  detector  is 
1.00  Mev. 

The  "averaged"  geometric  functions  of  these  detectors  were  obtained  by 
multiplying  the  transmission  functions  of  the  detectors  by  the  maximum  value 
of  the  geometric  function  of  the  detector  in  Channel  2  given  in  Figure  3-2h. 

The  results  are  depicted  in  Figure  ^-2b  and  listed  in  Table  Arl. 

3.3.7  COMPARISON  OF  CALIBRATION  RESULTS  WITH  INITIAL  ANALYSIS  OF  DETECTORS 

In  order  to  compare  the  experimental  values  of  the  geometric  functions 
of  the  detectors  with  those  which  were  computed  (SWC  50-TS-112)  for  the  pre¬ 
liminary  investigation  of  the  data,  the  geometric  functions  were  represented 
as  a  product  of  the  transmission  function  of  the  tube  and  the  geometric  tube 
factor.  The  function  was  obtained  experimentally  at  each  energy  E  by  dividing 
the  counting  rate  of  the  detector  by  the  value  of  the  planar  flux  of  electrons 
of  that  energy.  Hence,  the  geometric  factor  of  the  counter  was  taken  to  be 
the  geometric  function  for  E  >,>  E^,  where  E^  represents  the  threshold 
energy  of  the  counter.  The  transmission  function  of  the  tube  was  then  obtained 
by  dividing  the  geometric  function  by  the  geometric  factor.  The  following 
comparison  was  observed  between  the  experimental  values  for  Yj  «  90°,  a 

0°  and  the  values  computed  for  "Orientation  B"  given  in  SWC  50-TS-112. 

Channels  1  and  3?  The  measured  geometric  factor  was  about  30$ 
smaller  than  the  computed  value.  The  slope  of  the  transmission 
function  was  less  than  the  SWC-report  value  by  nearly  a  factor 
of  two.  (The  2g/cm  copper  plugs  in  windows  1,  5  and  6  diminished 
the  geometric  factor  by  a  factor  of  O.652  Instead  of  0. 5) • 

Channels  3  and  J:  The  geometric  factors  were  smaller  than  those 
given  in  the  SWC  report  by  18$.  The  transmission  functions  were 
in  fairly  good  agreement. 
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Channels,  2,  4  and  0:  The  measured  geometric  factors  were 
found  to  be  smaller  than  those  given  in  the  SWC  report  by 
about  4-0$.  The  transmission  functions  agreed  fairly  well. 

Channel  6:  The  empirical  method  which  was  used  to  determine 
the  calibrations  of  Channel  6,  described  in  Section 
is  considered  to  be  fairly  good  because  of  the  agreement 
obtained  by  testing  the  method  with  the  results  of  Channel  3* 

This  method  yielded  a  response  function  for  Channel  6  which 
is  much  lower  than  the  function  given  in  the  SWC  report  at 
energies  near  the  threshold  value,  being  smaller  by  a  factor 
of  ten  at  5  Mev  and  by  a  factor  of  two  at  10  Mev. 

Moreover,  the  initial  analysis  indicated  a  factor  of  tb»e  discrepancy 
between  the  flux  given  by  the  geometric  factors  of  the  long  counters  and 
that  given  by  the  short  detectors.  This  discrepancy  does  not  occur  when  the 
measured  geometric  functions  are  used. 
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ANALYSIS  OF  RESULTS 
4.1  PRESENTATION  OF  "RAW"  DATA 

The  "raw"  data  here  refer  to  the  counting  rates  of  the  detectors  as 
a  function  of  time.  The  computer  read-out  records  listed,  for  each  smple 
of  the  output  of  a  detector,  the  time  of  the  sample,  the  corrected  count¬ 
ing  rate,  and  the  orientation  angles  OJ  and  •  Of  these  counting  rates, 

only  those  which  were  sampled  when  occurred  in  the  interval  -20°^  Af 
20°  for  the  long  detectors  were  used;  and,  for  the  short  detectors  only  the 
counting  rates  which  were  sampled  when  the  collimator  axes  were  within 
+10°  of  the  iadiation  plane  were  used  in  the  analysis.  Data  which  were 
obtained  at  larger  displacements  were  discarded  because  the  geometric 
factors  of  the  detectors  began  to  change  rapidly  with  angle,  and  large  errors 
probably  would  have  ensued  due  to  the  non-planar  nature  of  the  actual  flux. 


Before  experimental  information  on  the  calibration  of  the  detectors 
was  obtained,  it  was  hoped  that  the  energy  dependence  of  the  geometric 
functions  of  the  long  detectors  would  vary  with  Yj  and  Aj^  so  that  a  more 
detailed  spectrum  could  have  been  obtained.  However,  essentially  no  change 
in  the  shape  of  the  geometric  functions  was  observed  in  the  angular  ranges 
69° <  ^  1110  and  aG  pointed  out  in  Section  3.  This  - 

orientation  angle  for  the  long  detectors  was  therefore  ignored. 


The  samples  of  the  counting  rates  which  occurred  within  the  angular  in¬ 
tervals  discussed  above  were  then  averaged  over  time  durations  for  which 
continuous  data  of  counts  per  pedestal  were  available  (5  to  12  seconds). 

The  process  was  straightforward  for  the  long  detectors:  the  counting  rates 
which  were  read-out  during  ifi  ^  20°  were  merely  added  and  divided  by  the 
total  number  of  samples.  However,  for  the  short  detectors,  it  was  necessary 
to  vary  the  Y  -angle  interval  over  which  the  counting  rates  were  accumulated 
in  a  mainer  which  was  dependent  on  Yj  .  This  was  accomplished  by  plotting 

two  curves  of  O^j  versus  ^  which  established  the  limits  of  these  angles 

for  a  collimator  displacement  from  the  radiation  plane  within  +10°.  The 
displacement,  S  ,  of  the  collimator  from  the  plane  is  given  by  the  equation, 

sin  cT  =  cos  60°  sin  "r|  Q  sin  -  jJ  )  -60°  ctn  v]  ■ 
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The  counting  rates  of  the  detectors,  averaged  as  discussed  above,  for 
the  five  rockets  which  passed  through  the  electron  shell  are  listed  in 


Tables  2  through  6  of  the  appendix, 
of  the  launch  times  of  the  rockets, 
times  and  the  burst  time  for  each  of 

Rocket  Flight 

2019 

2021 

202k 

2027 

2042 


The  tables  are  numbered  in  the  order 
The  time  interval  between  the  launch 
>  rockets  is  given  below: 

Time  of  Launch  after  Burst 

Hr: Min 

0:28 

1:58 
4:01 
18:42 
88:43 


After  averaging  some  of  the  counting  rates  given  in  the  tables  over  larger 
time  intervals  where  necessary  to  improve  the  statistical  accuracy,  the 
resulting  counting  rates  were  plotted  as  a  function  of  the  flight  time  of 
the  rockets.  These  data  are  shown  in  Figures  4-1  through  4-5,  again  in  the 
order  of  the  launch  times  of  the  rockets. 

The  errors  given  both  in  the  tables  and  in  the  figures  are  standard 
deviations  based  only  on  the  total  number  of  counts  used  in  the  averaging 
process.  At  high  counting  rates  the  major  errors  are  due  to  uncertainties 
in  the  determination  of  the  counts  per  pedestal  and  the  uncertainty  of  the 
"dead  time"  correction.  Systematic  errors  as  high  as  30%  or  more  may  have 
been  incurred  at  counting  rates  above  about  10^/sec.  The  probable  errors 
PR(J,  of  the  corrected  counting  rates  J^,  due  to  these  causes  is, 


%  ■  ^p/  +  Bo2pr' V<ltR0r). 

where 

R  m 

O 

observed  counting  rate,  15.75  N/Tc.  N  is  the  number 
of  counts  sampled  during  the  time  Tc/l5.75.  Tc  is 

the  period  of  the  commutator. 

PN  " 

probable  fractional  error  of  N  due  to  reading  of 

oscillograms. 

pr  - 

probable  error  of  dead  time  . 
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Fig.  4-1  Counting  Rates  of  Detectors  on  Flight  2019  as  Function  of  Flight  Time 
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Fig.  4-2  Counting  Rates  of  Detectors  on  Flight  2021  as  Function  of  Flight  Time 
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-3  Counting  Rates  of  Detectors  on  Flight  2024  as  Function 
of  Flight  Time 
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Fig.  4-4  Counting  Rates  of  Detectors  on  Flight  2027  as  Function  of  Flight  Time 
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Fig.  4-5  Counting  Rates  of  Detectors  on  Flight  2042  as  Function  of  Flight  Time 
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These  errors  are  listed  in  Table  4-1.  was  taken  to  be  10  microseconds. 


Table  4-1 

Probable  Errors  Due  to  Data  Readout  and 
Uncertainty  of  Dead  Time  Correction 


N 

( Counts ) 

PN 

RT 

(Counts/sec) 

pRj 

(Counts/i 

10 

0 

960 

10 

20 

0 

2060 

42 

30 

0 

3330 

230 

40 

.025 

4810 

4oo 

50 

.050 

6570 

64o 

60 

.075 

3680 

1560 

70 

.100 

11,300 

2380 

80 

.125 

14, 500 

4850 

90 

.150 

18, 700 

7350 

95 

.162 

21,300 

9180 

100 

.175 

24,300 

l4,4oo 

105 

.187 

27,900 

18,200 

110 

.200 

32,200 

23,400 

115 

.225 

37, 500 

30,600 

117 

.250 

40,000 

31*,1**) 

These  counting  rates  can  be  correlated  with  the  position  of  the  rockets 
in  space  as  well  as  with  the  locations  of  the  rockets  with  respect  to  geomag¬ 
netic  lines  of  force  with  the  information  given  in  Tables  7  through  11  of  the 
appendix  and  in  the  Figures  4-6,  4-7,  and  4-8.  These  tables  list,  against 
the  "flight  time"  of  a  rocket,  the  geographical  coordinates  of  the  rocket; 
the  magnetic  field  intensity,  B;  the  logarithm  of  the  integral  Invariant,  I; 
the  altitude  of  the  rocket  transformed  along  constant  local  B  and  I  values  to 
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ALTITUDE  (KM) 


NORTH  LATITUDE  (DEG) 

Fig.  4-6  Trajectories  of  Javelin  I  and  Jason  Rockets  Projected  in  Meridian  Plane  at  75° 
West  Longitude.  The  projection  of  the  trajectories  to  this  plane  was  done  by 
following  the  altitude  variations  of  local,  but  constant,  B  (geomagnetic  field 
intensity)  and  I  (integral  Invariant)  values  with  longitude.  Magnetic  lines  of 
force  which  intersect  the  surface  of  the  earth  at  every  degree  of  latitude  are 
also  shown.  The  electron  shell  was  centered  on  the  magnetic  line  of  force 
which  intersected  the  earth's  surface  at  a  latitude  of  about  33.6®. 
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o  * 

the  geographic  meridian  at  75  West  longitude,  H';  and  the  conjugate  mirror- 
point  altitude  (southern  hemisphere)  of  electrons  which  mirror  at  the  location 
of  the  rocket.  (The  last  column  will  be  referred  to  in  a  later  discussion. ) 

The  counting  rates  of  the  detectors,  when  compared  in  this  manner,  show 
that  the  peak  counting  rates  on  flights  2019>  2021,  2024,  and  2027  occurred 
when  the  trajectories  of  the  rockets  reached  a  shell  containing  the  magnetic 
line  of  force  which  intersects  the  earth's  surface  at  the  point  33*6°  North 
latitude,  75°  West  longitude.  The  pronounced  rise  to  the  peak  generally  began 
when  the  rockets  reached  the  shell  containing  the  magnetio  line  of  force 
which  falls  at  the  surface  of  the  earth  at  about  34.5°  North  latitude  along 
the  same  meridian.  High  counting  rates  also  occurred  at  greater  distances 
to  the  North  of  the  electron  shell  (the  wings).  The  counting  rates  of  all  of 
these  detectors  on  different  flights  but  at  comparable  spatial  locations 
decrease  with  time. 


The  B  and  I  values  used  for  the  analysis  described  in  this  report  were 
obtained  from  the  following  reference:  Tables  of  Adiabatic  Invariants 
for  Geomagnetic  Field  1955*0,  AFSWC-TN-60-8,  Kirtland  Air  Force  Base, 
New  Mexico,  i960. 
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4.2  BACKGROUND  RADIATION 
4.2.1  NATURAL  RADIATION 

A  large  effort  was  directed  towards  the  evaluation  of  the  background 
due  to  the  natural  trapped  electrons  for  each  of  the  flights.  The  result  of 
this  effort  was  less  than  satisfactory.  The  problem  was  greatly  complicated 
because  the  rockets  traversed  spatial  regions  where  the  natural  radiation  is 
known  to  be  highly  irregular.  In  general,  the  variation  of  this  radiation  is 
such  that  it  decreases  with  increasing  longitude,  it  increases  rapidly  with 
increasing  altltuue  along  magnetic  lines  of  force,  and  it  increases  even  more 
rapidly  across  magnetic  lines  of  force  towards  higher  latitudes. 

The  Javelin  I  rocket  data  (see  footnote  2  on  page  l),  which  yielded  the 
spectrum  and  intensity,  as  well  ds  pitch  angle  distributions,  of  the  trapped 
electrons  in  spatial  regions  similar  to  those  covered  by  the  Jason  rockets, 
verified  the  complexity  of  the  background  problem.  This  rocket  was  launched 
from  Wallops  Islsnd,  Virginia,  on  July  7,  1959,  about  10  months  after  the 
Argus  tests.  Its  trajectory  is  shown  in  Figures  4-6,  4-7,  and  4-8,  together 
with  those  of  the  Jason  rockets.  By  multiplying  the  spectral  intensity  of 
the  flux  obtained  with  the  Javelin  rocket  by  the  geometric  functions  of  the 
Jason  detectors  and  integrating  over  energy  (see  Equation  (l-l),  the  counting 
rates  of  the  Jason  detectors  along  the  Javelin  trajectory  were  inferred.  For 
the  purposes  of  this  discussion,  the  Jason  detectors  for  which  the  counting 
rates  were  computed  in  this  manner  will  be  referred  to  as  though  the  detectors 
were  actually  on  the  Javelin  rocket. 

A  comparison  of  the  counting  rates  of  the  detectors  on  the  Jason  rockets 
with  those  on  the  Javelin  revealed  that  in  the  wings,  a  high  fraction  of  the 
counting  rates  of  the  detectors  in  Channels  1  and  2,  and  a  lower  fraction  in 
Channels  7  6nd  8  were  due  to  the  natural  radiation.  In  fact,  a  plot  of  the 
counting  rates  of  the  detectors  in  Channel  1  for  the  five  Jason  flights  and 
the  Javelin  flight  as  a  function  of  altitude  revealed  that  the  counting  rates 
of  all  of  the  detectors  agreed  to  within  10^  along  the  ascent  trajectories 
to  an  altitude  of  600  km.  This  good  agreement  with  the  Javelin  data  was 
probably  accidental  since  the  absolute  measurement  of  the  flux  is  not 
expected  to  be  better  than  about  +  2 0$.  At  higher  altitudes,  a  separation  of 
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the  counting  rates  was  observed,  with  the  counting  rate  of  the  detector  cn 
the  Javelin  rocket  exceeding  the  rates  of  the  detectors  in  the  Jason  rockets 
because  of  the  continuation  of  the  Javelin  rocket  up  the  magnetic  line  of 
force . 

Efforts  to  determine  the  background  of  the  Jason  Channel  1  detector 
were  not  entirely  successful.  Estimates  of  the  background  were  made  by 
drawing  iso-count  lines  connecting  positions  on  the  ascent  and  descent  paths 
of  the  Javelin  rocket,  where  the  counting  rates  for  Channel  1  were  identical, 
and  assigning  these  rates  to  the  corresponding  Jason  detectors  at  the  spatial 
locations  where  these  lines  intersect  at  the  Jason  rocket  trajectories.  Such 
interpolations  were  attempted  using  the  graph  of  altitude  versus  latitude 
given  in  Figure  4-6,  and  the  graph  of  the  trajectories  on  the  B,  log  I  plot 
of  Figure  4-8.  The  resulting  background  counting  rates  of  the  detectors, 
when  plotted  against  the  flight  time  of  the  rockets,  gave  curves  which 
disagreed  somewhat  with  the  shapes  of  the  total  counting  rates  of  the  Channel 
1  detectors.  Moreover,  over  some  regions  of  the  curves  for  flights  2027  and 
2042,  the  background  counting  rates  exceeded  the  total  counting  rates  of  the 
detectors. 

It  was  clear  that  a  more  precise  computation  of  the  variation  of  the 
natural  radiation  was  required.  Various  empirical  approaches  were  attempted, 
based  on  the  variations  of  the  mirror-point  altitudes  in  the  southern  hemis¬ 
phere  of  the  electrons  detected  by  the  rockets.  A  comparison  of  the  altitude 
of  the  rocket  and  the  value  listed  as  Hs  in  Tables  A-7  through  A-12  shows  that 
the  mirror-points  altitudes  of  the  electrons  in  the  southern  hemisphere  were 
much  lower  than  those  at  the  rocket,  and  this  lowering  effect  was  enhanced  as 
the  rockets  moved  eastward.  The  results  gave  background  curves  which  again 
appeared  to  be  qualitatively  correct,  but  they,  too,  exceeded  the  total  count¬ 
ing  rates  in  some  areas. 

It  was  finally  decided  that  the  detectors  in  Channels  1,  5  and  2  could 
not  be  used  to  investigate  the  artificially  injected  electrons  in  the  regions 
of  the  wing  where  the  natural  radiation  was  intense. 
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Moreover,  the  accuracies  of  the  results  from  the  detectors  in  Channels 
7  and  8  were  seriously  impaired.  The  background  counting  rates  of  the 
Channel  7  detectors  computed  from  the  Javelin  data  were  about  a  factor  of 
two  lower  than  the  total  counting  rates  of  these  detectors  in  the  wings  on 
flights  2027  and  2042.  The  total  Channel  7  counting  rates  probably  included 
contributions  due  to  bremsstrahlung  from  the  natural  radiation  as  discussed 
below. 

Hence,  the  background  corrections  in  the  wings  for  the  Channel  7  detectors 
were  made  by  assuming  that  the  counting  rate  of  the  Channel  7  detector  on 
flight  2042  had  reached  the  background  level  and  that  this  natural  background 
level  was  the  same  for  all  of  the  flights.  The  rocket  trajectories  on  a  B-log 
I  plot  (Figure  4-8)  were  used  to  determine  these  backgrounds:  At  each  point 
of  a  trajectory,  the  background  was  obtained  by  means  of  a  log  I-weighted 
average,  at  constant  B,  of  the  Channel  7  counting  rates  measured  on  the  ascent 
and  descent  trajectories  of  the  2042  rocket.  This  evaluation  is  expected  to 
be  poor  for  the  portions  of  the  rocket  trajectories  which  depart  appreciably 
from  the  path  of  flight  2042.  Hence,  the  Channel  7  backgrounds  for  flight  2021, 
which  has  a  westward  velocity  component  (see  Figure  4-7),  and  for  flight  2024 
beyond  apogee  (see  Figure  4-8)  are  considered  to  be  poorly  known.  The  counting 
rates  of  the  Channel  7  detectors  on  flights  2019,  2021,  2024  and  2027,  together 
with  the  background  curves  given  by  this  method,  are  shown  in  Figures  4-9 
through  4-12,  respectively,  and  they  are  listed  in  Table  4-2.  The  probable 
errors  of  the  background  counting  rates  for  flights  2019,  2027,  and  a  portion 
of  2024  were  taken  to  be  20$,  the  extent  of  the  agreement  with  the  Javlin  I 
data  along  the  ascent  trajectories.  The  errors  were  arbitrarily  taken  to  be 
+50$  for  flight  2021  and  for  the  portion  of  flight  2024  beyond  apogee.  Table  4-2 
also  lists  the  net  counting  rates  of  the  Channel  7  detectors.  Total  probable 
errors  are  Included. 

4.2.2  BACKGROUND  DUE  TO  PREVIOUSLY  DETONATED  DEVICES 

Comparison  of  the  counting  rates  of  the  Channel  1  detectors  on  Patrick 
launched  flights  1822,  1909,  2020,  and  2043  at  similar  positions  of  the 
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Fig.  4-9  Counting  Rate  of  Detector  in  Channel  7  on  Flight  2019  and  Background 
Correction  for  This  Detector 
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Fig.  4-10  Counting  Rate  of  Detector  in  Channel  7  on  Flight  2021  and  Background 
Correction  for  This  Detector 
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Fig.  4-11  Counting  Rate  of  Detector  in  Channel  7  on  Flight  2024  and  Back¬ 
ground  Correction  for  This  Detector 
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-12  Counting  Rate  of  Detector  in  Channel  7  on  Flight  2027  and  Background 
Correction  for  This  Detector 
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trajectories  of  these  rockets  on  a  B,  log  I  plot,  revealed  an  abnormally  high 
counting  rate  for  the  earliest  flight  1822.  This  high  counting  rate  may  have 
been  caused  by  the  presence  of  electrons  injected  into  the  geomagnetic  field 
by  the  detonation  of  nuclear  devices  at  high  altitudes  in  the  Pacific,  as 
was  pointed  out  in  the  initial  analysis  of  these  data.  The  decay  rate  of 
these  electrons  was  so  great,  however,  that  the  background  due  to  this  cause 
was  considered  to  be  negligible  for  the  later  flights  now  under  discussion. 

(For  the  few  rockets  launched  from  Patrick  and  Ramey  for  which  the  body  motions 
of  the  rockets  were  estimated,  the  counting  rates  of  the  small  detectors  were 
incorrectly  analyzed  in  the  preliminary  report  because  the  collimator  axes  of 
these  detectors  never  entered  the  central  plane  of  the  radiation.) 

4.2.3  COSMIC  RAY  BACKGROUND 

The  counting  rates  produced  in  the  detectors  by  cosmic  rays  were  evaluated 
from  the  counting  rates  of  the  tubes  which  had  high  threshold  energies  and  in 
regions  of  the  flights  where  the  bomb  and  natural  radiations  were  low.  Hence, 
the  background  counting  rate  of  the  long  tubes  was  taken  to  be  the  mean  count¬ 
ing  rate  of  the  Channel  6  detector  observed  on  flight  1822,  which  was  launched 
prior  to  a  detonation  of  an  Argus  device  and  into  the  slot  between  the  natural 
radiation  zones.  Thi3  counting  rate  was  found  to  be  16  +  2  counts/sec.  It 
also  agreed  with  the  mean  counting  rate  obtained  for  the  Channel  6  counting 

rate  on  flight  2027  and  with  the  rate  computed  by  using  the  omnidirectional 

2 

geometric  function  of  the  detector  and  a  cosmic  ray  flux  of  1.5  cm  sec. 

The  counting  rate  of  the  small  detectors  due  to  cosmic  ray3,  determined 
by  the  method  indicated  above,  also  agreed  with  the  background  counting  rates 
determined  1.  r  the  same  geiger  tube  types  (Anton  222R)  used  on  the  Javelin  1 
rocket.  This  counting  rate  was  found  to  be  2.1  +  0.25/sec. 

4.2.4  BREMSSTRAHLUfJG 

There  was  early  evidence  that  the  Channel  6  detector  was  responding  to 
bremsctrahlung  in  the  wing.  The  counting  rate  of  this  detector,  above  the 
cosmic  ray  background,  yielded  curves  when  plotted  against  flight  time  which 
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were  quite  different  in  shape  from  those  given  by  the  other  detectors  on  the 
same  flight.  The  Channel  6  curves  were  more  peaked  in  the  wing  where  the 
natural  radiation  was  more  intense.  Moreover,  the  decay  rate  of  the  outputs 
of  the  Channel  6  detectors  in  the  wing,  obtained  by  plotting  the  counting 
rates  of  these  detectors  above  cosmic  ray  background  as  a  function  of  the 
launchtime  of  the  rockets  after  the  burst,  was  found  to  be  much  smaller  than 
that  given  by  the  detectors  in  Channels  3  and  7« 

A  rough  calculation  was  performed  to  evaluate  the  response  of  this 
Channel  6  detector  to  the  bremsstrahlung  of  the  natural  radiation.  The 
detector  response  was  computed  for  the  electron  spectrum  given  by  the  Javelin 
I  measurements,  and  the  natural  radiation  level  at  the  position  of  a  particu¬ 
lar  Jason  rocket  was  estimated  by  integrating  over  energy  the  product  of  the 
Javelin  I  spectrum  and  the  geometric  function  of  the  Channel  1  detector.  It 
was  found  that  the  ratio  of  the  counting  rate  of  the  Channel  6  detector  due  to 
this  bremsstrahlung  to  the  counting  rate  of  the  Channel  1  detector  was  about 
0.001  +  0.0005.  The  maximum  correction  amounted  to  twice  cosmic  ray  back¬ 
ground. 

Section  4.3  NET  COUNTING  RATES  OF  CHANNEL  7,  3,  and  6  DETECTORS 

The  background  correction  for  Channel  7»  obtained  as  described  in 
Section  4.2.1,  and  the  net  counting  rates  of  the  detectors  in  Channels  7>  3, 
and  6  for  flights  2019,  2021,  2024,  and  2027  are  listed  in  Table  4-2.  For 
flight  2042,  the  true  counting  rate  of  the  Channel  7  detector,  which  was  used 
to  compute  background  corrections  for  the  other  flights,  and  the  net  counting 
rate  of  the  Channel  3  detector  are  listed  in  this  table. 

Only  the  cosmic  ray  background  counting  rate  was  deduced  from  the  true 
counting  rates  of  the  Channel  3  detectors.  The  bremsstrahlung  correction 
described  in  Section  4.2.4,  as  well  as  the  cosmic  ray  background  correction, 
were  applied  to  the  Channel  6  readings. 
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5.  RESULTS  AND  DISCUSSION 

Many  limitations  of  the  Jason  experiment,  some  of  which  were  discussed 
in  Sections  4.1  and  4.2,  make  it  difficult  to  draw  definitive  conclusions 
concerning  the  electron  injection  and  trapping  phenomena  which  occurred 
following  the  Argus  II  detonation.  However,  the  more  extensive  analysis 
done  in  the  course  of  the  present  contract  does  allow  critical  examination 
of  many  aspects  of  the  problem.  Wherever  possible,  definite  conclusions 
have  been  drawn  and  are  presented  in  this  section.  However,  in  many  cases 
the  data  are  insufficient  to  establish  the  existence  of  expected  effects. 

In  particular,  some  of  the  conclusions  reached  after  the  preliminary  analysis 
are  not  borne  out  by  the  present  work,  after  due  allowance  has  been  made 
for  experimental  uncertainties. 

The  major  difficulty  in  obtaining  a  complete  picture  of  the  electron 
trapping  phenomena  from  the  Jason  flights  results  directly  from  the  scar¬ 
city  and  sporatic  nature  of  experimental  data.  The  five  flights  which  have 
been  examined  in  detail  passed  through  the  band  at  different  times,  at 
different  altitudes,  and  at  different  longitudes.  Since  the  electron  fluxes 
are  functions  of  position  and  time  after  the  burst  and  since  the  geomagnetic 
field  is  not  regular  in  this  region,  it  is  not  clear  how  the  data  from 
different  flights  can  best  be  compared  to  obtain  precise  information,  A 
second  deficiency  of  the  data  results  from  the  restricted  dynamic  range  of 
the  detectors,  so  that  only  a  few  of  the  counters  gave  information  at  the 
same  time.  In  most  cases  where  the  small  counters  had  sufficient  counting 
rates  for  statistical  accuracy,  the  large  counters  were  counting  so  rapidly 
that  the  telemetering  tracings  could  not  be  read  or  the  uncertainties  in  the 
dead  time  corrections  were  excessive.  This  difficulty  severely  limits  a 
comparison  of  the  energy  spectrum  obtained  during  various  flights,  and  lowers 
the  precision  with  which  changes  in  the  spectrum  can  be  detected.  A  final 
serious  deficiency  was  in  the  background  counting  rate  information.  This 
point  has  been  discussed  in  Section  4.2.  However,  it  should  be  noted  here 
that  during  most  of  the  time  when  the  rockets  were  not  in  the  intense  portion 
of  the  radiation  band,  background  uncertainties  constituted  the  largest 
source  of  error. 
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5.1  GENERAL  RESULTS 

The  results  of  the  analysis  are  summarized  in  Figures  4-1  through  4-5, 
which  give  the  counting  rates  as  a  function  of  time  after  launch  of  flights 
2019,  2021,  2024,  2027  and  2042.  In  addition  to  the  localized  region  of  high 
IntaiBity  radiation  which  appears  at  the  far  right-hand  side  of  the  figures, 
there  is  a  broad  region  at  higher  latitudes  where  the  radiation  is  signifi¬ 
cantly  above  the  background  levels.  This  feature  is  most  pronounced  for 
Channels  3  and  6,  but  is  also  apparent  for  Channel  J.  Because  of  the  large 
background  of  low  energy  Van  Allen  electrons,  the  enhanced  radiation  in  this 
wing  i3  not  as  apparent  for  the  low  energy  channels.  The  difference  in  the 
magnitude  of  ..he  flux  in  the  intense  band  for  the  various  flights  is  due 
primarily  to  the  decay  with  time  of  the  shell  of  trapped  electrons,  although 
this  effect  is  modified  considerably  by  the  differences  in  the  altitudes  at 
which  the  various  vehicles  passed  through  the  band.  For  example,  flight  2024 
which  passed  through  the  band  at  an  altitude  of  400  km  showed  a  more  pronounced 
peak  in  counting  rate  than  did  flight  2021,  which  actually  took  place  earlier, 
but  crossed  the  band  at  an  altitude  of  only  about  100  km.  These  altitudes 
refer  to  the  heights  of  the  conjugate  points  of  the  rocket  positions. 

5.2  VARIATION  OF  ELECTRON  ENERGY  SPECTRUM  WITH  POSITION 

The  ratios  of  the  counting  rates  of  the  detectors  were  examined  at 
various  times  during  each  flight  to  see  if  the  spectrum  of  the  injected  electrons 
was  different  at  different  points  of  the  trajectory.  In  particular,  it  was  of 
interest  to  see  if  the  spectrum  in  the  wing  differed  from  that  of  the  intense 
band,  since  this  would  shed  light  on  whether  the  wing  electrons  were  injected 
directly  into  their  orbits  or  whether  they  diffused  across  field  lines.  From 
the  discussion  of  the  background  problem  in  Section  4.2,  it  is  aoparent  that 
only  the  counting-rate  ratios  R^/R^  and  R_/R-  can  be  considered  reliable  for 
this  purpose.  These  ratios  for  flights  2019  and  2024  are  listed  in  Table  5-1, 

In  the  center  of  the  band,  both  detectors  in  Channels  7  and  3  on  flight  2019 
were  saturated.  The  value  listed  in  parenthesis  for  this  flight  in  the  \j/ R 
column  is  the  9,/R,  ratio  which  was  measured  in  the  center  of  the  band.  The 
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listing  in  parenthesis  in  the  R^/Rg  column  was  obtained  by  integrating  the 
product  of  the  geometric  function  of  the  Channel  3  detector  and  the  best  fit 
spectrum  over  energy  and  dividing  the  result  by  the  measured  rate  Rg.  The 
ratios  in  the  table  appear  to  be  changing  systematically  with  time,  with  the 
R^/R^  indicating  a  softening  of  the  spectrum  and  the  R^/Rg  ra^°  indicating 
a  hardening  of  the  spectrum.  It  is  not  believed  that  these  ratios  demonstrate 
a  spectrum  change.  It  is  more  likely  that  the  trend  reflects  a  systematic 
error  in  the  evaluation  of  the  background. 


TABLE  $-1 

Comparison  of  Counting-Rate  Ratios 


Flight  Time 
(sec) 

(r7t  -  Kfm*  -  r,b) 

t«3T  -  r3b>/cr6T  •  »6B> 

Flight  2019  Flight  2024 

Flight  2019  Flight  2024 

250 

I.56  +  .26  1.60  +  .80 

300 

1.46  +  .28  1.81  +  .74 

350 

1.50  +  .33  1.86  +  .86 

720  +  1600 

1*00 

1.50  +  .35  1.75  ±  .92 

660  +  1430 

450 

1.60  +  .40  2.14  +  .82 

520  +  900 

500 

1.70  +  .45  2.77  +  .86 

390  +  500 

550 

1.86  +  .49  3.50  +  1.50 

310  +  310 

600 

1.92  +; .53  3.00  +  .80 

290  +  290 

650 

2.04  +  .57  1.85  +  .71 

250  +  200 

700 

2.25  +  .72 

240  +  150 

Band 

(1.90  +  .47)* 

(4l6  +  136)*  245  +  132 

* 


Refer  to  text  for  comments  on  this  ratio. 
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5.3  VARIATION  OF  ELECTRON  ENERGY  SPECTRUM  WITH  TIME 

The  reduced  data  were  also  examined  carefully  to  see  if  there  was  an 
indication  of  a  change  in  the  electron  energy  distribution  during  the  88- 
hour  period  covered  by  the  Jason  experiment.  Since  no  definitie  variation 
of  spectrum  with  position  had  been  detected,  it  was  felt  possible  to  compare 
data  taken  on  different  flights  without  regard  to  the  position  of  the  vehicle 
during  the  time  the  comparison  was  made.  In  this  way  sections  of  various 
flights  could  be  compared  in  which  the  counting  rates  were  similar  and  dead¬ 
time  corrections  would  be  comparable.  This  investigation  failed  to  show  any 
definite  time  variation  in  the  spectral  distribution,  in  agreement  with  the 
preliminary  analysis  presented  two  years  ago.  Unfortunately,  the  low  energy 
channels  where  the  spectrum  change  is  expected  to  be  most  pronounced  were  not 
very  useful  in  this  regard  because  of  the  high  and  uncertain  background  pro¬ 
duced  by  the  Van  Allen  electrons.  The  fact  that  the  spectrum  obtained  during 
flight  2019  which  took  place  about  half  an  hour  after  the  explosion  was  not 
significantly  different  from  the  later  flights  is  perhaps  an  indication  of  the 
poor  sensitivity  of  the  spectrum  test. 

5.4  DECAY  OF  ELECTRON  FLUX 

A  quantitative  measurement  of  the  decay  of  geomagnetically  trapped 
electrons  was  one  of  the  most  Important  objectives  of  the  Argus  experiment. 

The  observed  decay  time  of  several  days  is  considerably  less  than  the  decay 
time  of  many  years  which  has  been  calculated  for  particles  distributed  in  the 
Van  Allen  electron  belt.  The  relative  short  time  of  decay  has  been  attributed 
to  the  low  injection  altitude  and  to  the  low  altitude  at  which  the  measurements 
were  made.  It  has  been  argued  that  what  is  observed  is  the  decay  of  a  higher 
order  normal  mode  in  the  spatial  distribution  of  electrons  in  the  trapping 
volume;  hence  the  lifetime  is  correspondingly  reduced.  Theoretical  calcula¬ 
tions  have  suggested  that  the  intensity  of  this  distribution  may  decrease  as 
the  inverse  first  power  of  the  time  after  detonation,  a  time  variation  which 

is  in  fair  agreement  with  observation.  However,  it  has  also  been  observed  that 

-1.2 

the  beta  decay  of  gross  fission  fragments  behaves  as  t  .  Since  the 
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interpretation  which  would  be  placed  on  the  Argus  experiment  is  vastly  differ¬ 
ent  depending  on  which  time  behavior  is  obtained,  a  considerable  effort  was 
made  to  resolve  this  question  from  the  results  of  the  Jason  rockets. 

Because  of  the  different  trajectories  followed  by  the  various  flights 
and  because  no  method  of  unquestioned  validity  was  found  to  extrapolate  the 
measurement  from  one  position  to  another,  it  was  not  possible  to  compare  in¬ 
tensities  found  at  different  times  near  the  center  of  the  band.  However,  in 
the  broad  region  of  radiation  north  of  the  band,  it  is  possible  to  make  a 
much  more  accurate  estimate  of  the  decay  rate.  From  the  curves  of  Figure  4-8 
which  give  a  plot  of  the  vehicle  trajectories  on  B,  Log  I  coordinates,  it  can 
be  seen  that  there  are  seven  points  where  two  of  the  trajectories  Intersect. 

At  these  points  the  rockets  of  the  two  flights  pass  through  equivalent  points 
on  the  magnetic  shell  and  a  comparison  of  the  counting  rates  of  the  two  flights 
will  then  give  an  estimate  of  the  flux  decay  during  the  time  Interval  between 
the  two  flights.  A  small  correction  was  made  for  the  difference  in  the  longi¬ 
tudes  of  the  two  vehicles  at  the  time  of  the  intersection.  This  correction 
is,  in  essence,  the  variation  of  flux  along  a  line  of  constant  B  and  Log  I. 

It  was  obtained  from  flight  2025  launched  from  Cape  Canaveral.  This  flight 
was  directed  towards  the  east,  and  the  ascending  and  descending  portions  of 
the  trajectory  passed  through  equivalent  magnetic  locations,  which  were  dis¬ 
placed  in  longitude  and  also  in  altitude.  Hence,  this  flight  gives  the  rate 
at  which  the  flux  varies  during  the  eastward  drift.  Corrections  for  this 
effect  sure  fairly  small,  ranging  from  0  to  7  percent  because  of  the  high 
altitude  at  which  most  of  the  intersections  occurred. 

Since  each  intersection  gives  a  value  of  the  decay  during  a  particular  time 
Interval,  the  seven  Intersections  give  seven  points  on  a  time  decay  curve. 

These  points,  obtained  for  the  Channel  3  detector  are  shown  in  Figure  5-1. 

Here,  the  ratio  of  the  counting  rates  observed  by  the  two  rockets  at  each  inter¬ 
section  point  is  plotted  against  the  ratio  of  the  times  at  which  the  rockets 
reached  the  point  of  intersection  after  the  time  of  burst.  The  error  bars 
indicate  the  total  probable  error  which  results  from  counting  statistics, 
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dead-time  uncertainties,  background  uncertainties,  and  trajectory  uncertain¬ 
ties.  Assuming  that  the  decay  has  a  function  form  of  t  n,  a  least  square  fit 
to  the  data  gives  a  value  of  n  of  1.15  +  .03*  Since  the  ratio  of  the  counting 
rates  must  be  equal  to  1  when  the  time  ratio  equals  1,  only  a  one  parameter 
fit  (the  slope  of  the  line)  was  required.  The  same  calculation  for  channel  7 
gives  similar  results  but  with  considerably  greater  probable  errors. 

On  the  basis  of  this  analysis,  it  appears  more  likely  that  the  time 
dependence  follows  that  of  gross  fission  product  decay,  rather  than  the 
predicted  t”b  The  line  of  slope  n  =■  1  is  drawn  in  Figure  4-13  for  comparison. 

5.5  ENERGY  SPECTRUM  OF  ELECTRONS 

Since  no  variation  in  time  or  space  was  noted  in  the  energy  distribution 
of  the  electrons  in  the  band,  the  determination  of  the  electron  spectrum  was 
made  at  positions  where  the  experimental  data  were  most  reliable.  The  best 
Information  for  this  purpose  was  obtained  in  flight  2019  during  its  passage 
through  the  band,  at  which  time  reliable  data  were  obtained  simultaneously 
with  detectors  2,  4,  6  and  8.  The  spectrum  was  represented  by  the  expression 

f(E)  =  AE0  exp(-BE) 

where  E  is  the  electron  energy  in  Mev  and  A,  B,  and  D  are  parameters.  Since 
A  merely  normalizes  the  flux,  the  shape  parameters  D  and  B  are  of  greatest 
concern. 

To  deduce  the  electron  spectrum,  a  set  of  values  of  B  and  D  were  selected 
and  the  counting  rates  calculated  from  the  equation, 

Ri  =  fG i(E)  f  (E)  dE 

where  G^E)  is  the  geometric  factor  of  che  "i"th  detector  as  obtained  by  the 
calibration  measurements  described  in  Section  3.  The  calculated  counting  rates 
were  then  compared  to  the  experimental  values  by  forming  the  sum  of  the  squares 
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of  the  deviations 

M(B,D)  -  \A/±  (R£  -  Ri«x)2 

where  is  a  weighting  function  equal  to  the  reciprocal  of  the  square 

of  the  probable  error  in  the  counting  rates,  and  R^ex  is  the  experimental 
counting  rate  of  the  "i"th  detector.  The  values  of  M  form  a  matrix  in  B,  D 
space  and  the  minimum  value  of  M  as  a  function  of  B  and  D  indipates  the  values 
of  B  and  D  which  best  represent  the  experimental  result*  This  computing 
process  was  carried  out  using  a  digital  computer.  The  resulting  spectrum 
was  modified  by  interation  to  improve  the  agreement  with  the  counting  rates. 

The  best  fit  spectrum  is  depicted  by  the  solid  line  in  Figure  5-2.  An  indi¬ 
cation  of  the  accuracy  of  the  spectrum  is  given  by  the  rectangles  shown  in 
the  figure.  These  rectangles  are  the  "raw -data"  defined  as  follows:  the 
vertical  lines  of  the  rectangles  are  drawn  at  the  energy  values  which  are 
at  the  half-height  points  of  the  product  of  the  best  fit  spectrum  and  the 
geometric  functions  of  the  detectors.  The  vertical  positions  of  the  centers 
of  the  rectangles  are  located  at  the  geometric-function  weighted  averages  of 
the  flux  in  these  energy  intervals;  and  the  heights  of  the  rectangles  corres¬ 
pond  to  twice  the  probable  errors  of  the  net  counting  rates  of  the  detectors. 

Two  other  spectra  are  shown  in  Figure  $42  for  comparison:  An  unmodified 

fission  spectrum  and  a  fission  spectrum  divided  by  the  product  CO,  T_,  the 

d  b 

azimuthal  drift  velocity  times  the  bounce  period  between  mirror  points.  The 
disagreement  of  the  best  fit  spectrum  with  the  fission  spectrum  is  not  severe; 
it  is  now  much  more  in  agreement  with  the  data  than  was  indicated  by  the 
Initial  analysis  )f  the  data.  However,  the  better  agreement  of  the  data  with 
the  fission  spectrum  divided  by  CO^  T^  might  be  evidence  of  a  fission  spectrum 
source,  a  lifetime  of  less  than  a  drift  period  for  the  electrons  which  reach 
the  region  of  the  Jason  rockets,  and  no  evidence  of  energy  degradation  by 
scattering. 

This  statement  can  be  verified  as  follows:  the  omnidirectional  flux  given 
by  the  Jason  rockets  is  equal  to  the  number  density  of  the  electrons  times 
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Fig.  5-2  Omnidirectional  Flux  at  Center  of  Electron  Shell 
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their  velocity.  Or, 


f(E)dE  ■  n(E)dE  .  v,  (5-l) 

where  n(E)  is  the  number  density  of  electrons  of  energy  E  per  unit  energy 
interval.  If  N(E)  designates  the  number  of  electrons  emitted  by  the  fission 
spectrum  per  unit  time  and  unit  energy  interval  and  dV  is  the  incremental 
volume  the  electrons  fill  during  the  incremental  emission  time  dt,  then 

f  ( E)dE  -  -^V-E-d-  •  v  (5-2) 

This  equation  shows  that  a  fission  spectrum  multiplied  by  electron  velocity 
will  be  measured  if  there  is  no  electron  loss  and  if  the  electrons  are 
injected  into  a  magnetic  shell  of  fixed  volume  for  a  period  of  time  which  is 
long  in  comparison  to  drift  periods.  However,  if  the  electrons  which  were 
reflected  in  the  region  of  the  Jason  rockets  become  lost  beyond  the  longitude 
of  the  measurements,  then  for  these  electrons  dV/dt  is  proportional  to 
and  v  is  inversely  proportional  to  T^.  Hence,  with  no  loss  of  energy  by 
scattering,  the  measurements  should  yield  a  fission  spectrum  divided  by 

^d  Tb* 

The  observed  counting  rates  R.  of  the  detectors  in  the  center  of  the  band 

F  * 

on  flight  2042  and  the  rates  R^  given  by  the  product  of  the  modlfled-fission 
beta  spectrum  ( N(  E)/ ^  T  )  and  the  geometric  functions  of  the  detectors  inte- 


grated  over  energy  are  listed  below 

Channel  Number  2 

for  comparison. 

a 

4 

6 

1+030  +  760 

2000  +  360 

1050  +  190 

1280  +  2)0 

4030 

2730 

1000 

1250 

For  flight  2019,  which  was  launched  about  28  minutes  after  the  burst,  the 
good  agreement  of  the  measurements  with  this  spectrum,  is  exceedingly  difficult 
to  understand.  The  mechanism  proposed  earlier  involved  the  passage  of  fission 
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fragments  from  the  detonation  point  to  the  conjugate  region  via  the  magnetic 
field  lines,  and  electrons  which  were  emitted  during  the  time  the  fission 
fragments  were  in  transit  would  be  trapped.  If  this  process  did  occur,  the 
injection  would  take  place  only  during  the  fragment  passage  and  the  spectrum 
observed  during  flight  2019  would  consist  of  a  line  spectrum  containing  con¬ 
tributions  of  particles  whose  energies  were  appropriate  for  drift  around  the 
world  one  or  more  times  during  the  interval  between  the  detonation  time  and 
the  measurement.  Comparison  of  the  counting  rates  <i£)  expected  with  this 
spectrum  with  the  observed  rates  R^,  as  shown  in  Table  5-2,  reveals  that  the 
agreement  is  poor. 

An  alternate  hypothesis  for  injection  would  be  to  assume  that  the  frag¬ 
ments  become  distributed  in  such  a  manner  that  injection  can  take  place  for 
a  much  longer  period  of  time.  If  this  were  the  case,  the  spectrum  observed 
in  flight  2019  would  be  continuous,  containing  components  for  very  high 
energies  but  being  cut-off  at  energies  below  about  0.8  Mev,  since  electrons 
below  this  energy  cannot  drift  the  required  distance  in  the  available  time. 

The  observed  counting  rates  (R^)  and  the  rates  (R^)  calculated  with  these 
assumptions  are  also  listed  for  comparison  in  Table  5-2.  In  this  case,  the 

high  threshold  energy  detector  in  Channel  6  reads  quite  low  because  of  the 

•1.2 

t  *  decay  of  the  source  (these  high-energy  electrons  were  emitted  at  a  later 
time  than  the  low-energy  electrons  to  have  arrived  at  the  rocket  position  at 
the  seme  time). 

It  should  be  noted  that  the  above  calculations  assume  that  no  energy  is 
lost  by  the  electrons  during  their  drift  about  the  earth.  It  is  unlikely 
that  the  electrons  do  not  loose  energy  since  a  considerable  amount  of  scatter¬ 
ing  must  occur  if  the  electrons  are  to  be  found  in  orbits  which  are  accessible 
to  the  Jason  rockets.  Due  to  the  anamolies  in  the  geomagnetic  field,  electrons 
whose  mirroring  points  are  above  the  earth's  surface  on  the  magnetic  line  pass¬ 
ing  through  the  detonation  point  will,  after  drifting  to  the  latitude  of  the 
Jason  flights,  mirror  at  altitudes  far  above  apogee  of  the  vehicles.  Therefore, 
in  order  to  be  detected  these  electrons  must  have  been  scattered  into  lower 
altitude  orbits,  and  their  energy  must  be  correspondingly  modified.  The  amount 
of  energy  loss  expected  cannot  be  calculated  under  the  present  state  of  the 
trapping  theory  and  the  uncertainty  introduced  by  the  possibility  of  energy 
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loss  makes  the  conclusions  derived  from  spectrum  information  extremely  tenta¬ 
tive. 


TABLE  5-2 


Comparison  of  Counting-Rate  Ratios  Based  on  the 
Various  Injection  Mechanisms  for  Flight  2019 


CHANNEL  NUMBER(i) 

2  8 

k 

6 

Observed  Ratios  1.00 

(Ri/R2) 

+  .27  0.50  +  .13 

0.2 6  +  .07 

0.32  +  .09 

Modified  Fission 

Beta  Spectrum  (R^/R^) 

1.00  0.68 

0.25 

0.31 

Line  Spectrum  (R^/R^) 
(instantaneous  injection) 

1.00  0.88 

0.18 

0.19 

Continuous  Ehtission 

1.00  0.90 

0.22 

0.042 

(rC/rC)  (Cut-off  below  .8 
Mev  because  of  drift  line 
requirements) 
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5.6  INJECTION  EFFICIENCY 

A  lower  limit  of  the  trapping  efficiency  can  be  determined  by  consid¬ 
ering  the  Jason  data  to  be  a  measure  of  the  rate  at  which  electrons  are  lost. 
This  is  plausible  because  of  the  special  region  in  which  the  data  were  taken 
with  respect  to  the  total  configuration  of  the  geomagnetic  field.  As  discussed 
earlier,  electrons  which  reached  the  region  of  the  Jason  flights  were  very 
likely  absorbed  by  the  atmosphere  as  they  proceeded  to  drift  eastwards.  Hence, 
an  integration  of  the  •'’lux  measurements  over  time  gives  a  lower  limit  of  the 
total  loss  c.f  trapped  electrons,  from  which  a  lower  limit  of  the  trapping 
efficiency  can  be  estimated. 

The  incremental  loss  rate  of  electrons,  observed  by  the  Jason  rockets, 
can  be  obtained  by  finding  the  number  of  electrons  in  an  incremental  magnetic 
6hell  which  cross  a  certain  area  at  the  equator  per  unit  time  and  which  have 
pitch  angles  at  the  equator  such  that  the  electrons  reach  the  altitudes  of 
the  Jason  rockets.  By  Liouville's  theorem,  we  know  that  the  electron  flux  at 
the  equator,  directed  at  the  angle  0 (  with  respect  to  the  geomagnetic  line 

of  force,  is  equal  to  the  flux  at  the  Jason-rocket  altitudes  which  is  directed 
essentially  at  90°  to  the  magnetic  field  line.  The  angle  Oi  can  be  obtained 

®  O 

from  the  adiabatic  invariance  of  the  magnetic  moment  of  the  electrons  (sin' OC 

o 

equals  ratio  of  magnetic  field  intensity  at  equator  to  magnetic  field  intensity 
at  Jason  altitude  along  the  same  magnetic  line  of  force).  Thus,  the  number  of 
electrons  of  energy  E  in  the  range  E  which  cross  a  unit  cross  sectional  area 
of  the  incremental  magnetic  shell  at  the  equator  at  the  time  t  is  f(E,  t)  /\ E/ 
cos  0(  q,  where  f(E,  t)  is  the  measured  flux. 

To  get  the  number  of  electrons  which  are  lost  per  second  from  the  shell, 
the  flux  f(E,  t)  ^  E/cos  must  be  multiplied  by  the  cross  sectional  area 
rQ  ArQ  A  of  ^e  magnetic  shell.  Here  rQ  is  the  geocentric  distance 
to  the  line  of  force  at  the  magnetic  equator,  £  rQ  is  the  thickness  of  the 
shell,  and  ^  ^(E)  *»  T^/2.  Note  that  if  the  azimuthal  extent  h*  of 

this  "loss  area"  were  taken  to  be  larger  than  CO^  T^/2,  then  electrons  would 

be  ’fcounted"  more  than  once;  if  A  ^  were  taken  to  be  smaller  than  CO.  T  /2. 

1  d  b'  9 

then  some  of  these  electrons  could  avoid  striking  the  "loss  area"  if  their 


103 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  ond  SPACE  DIVISION 


2-96-61-1 


phases  were  appropriate*  Hence,  the  loss  of  electrons  of  energy  E  in  the 
range  &  E,  contained  in  the  shell  at  rQ  of  thickness  £  rQ  per  second  is 


L(E,ro,t)  AE  ^rQ 


f(E,t)  Ae 

COB  Q(  _ 

o 


r0  Aro  (5-3) 


In  the  previous  section  it  was  shown  that  the  measured  spectrum  agreed 
exceedingly  well  with  a  fission  spectrum  divided  by  Ck)^  T  •  Consequently, 
since  the  measured  spectrum  in  Equation  (5-3)  1®  multiplied  by  At^®)  "  ^d^b^2’ 
the  loss  spectrum  L(E)  is  Just  the  fission  spectrum. 

An  evaluation  of  L(E)  in  the  center  of  the  band  at  t  »  4l  minutes  (2019 
data)  and  a  subsequent  integration  over  energy  yields  the  following  loss  rate 
of  electrons  from  the  band: 

flO  Mev 

fc0(4lmin.)  =  L(E),  2re,  4l)  (.047  re)  dE 

/.  5  Mev 

lfl 

»  9.44  x  10J'0  electrons/sec. 


A  similar  evaluation  and  integration  of  Equation  (5-3)  over  the  "north 
wing"  of  the  trapping  region  yields, 

a  i  O 

^NW  ^3^  ^  t  ^  4o.7  min)  a  4.23  *  10  electrons/sec. 

Now,  by  assuming  that  the  loss  rate  from  the  low  latitude  wing  is  equal  to 

that  from  the  high  latitude  wing  and  extrapolating  the  time  decay  of  the  flux 

(  t“l*15j  which  was  observed  in  the  range  32  min  ^  t  ^  5333  min,  to,  say, 

t  a  1  min,  the  total  loss  of  electrons  after  t  a  l  min,  becomes, 

-  eo 

N°  ^  Jl  tV4l)  +  2Nnw(4i)](i^)  dt>  or 

21 

N  6.2  x  10  electrons, 

o 

Finally  by  taking  the  toted  number  of  fission  betas  emitted  one  minute  or  later 
after  the  detonation  to  be 
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ML  -  fl.7  KT)  (1.3  x  1023  fiBa1^  frag>-)  (2.1  electrons/fission 
T  emitted  after  one  minute) 

■  4.65  x  10  J  electrons, 
the  lower  limit  of  the  trapping  efficiency  is, 


(Trap.  Rff.)lwer  llalt  -  6.2  x  1021A.65  x  1023  -  1.3  x  10-2  . 


This  value  is  believed  to  be  a  lover  limltbecause  it  is  quite  likely  that 
more  electrons  were  lost  than  those  detected  by  ttye  Jason  rockets.  Also,  this 
trapping  efficiency  is  on  the  low  side  because  no  account  was  taken  of 
electrons  which  might  have  been  trapped  for  very  long  periods  of  time,  that  is, 
of  electrons  with  mirror  points  so  high  as  to  remain  trapped  during  the  course 
of  this  experiment. 


5.7  COMMENTS  ON  RESULTS 

The  observed  decay  rate  of  the  electron  flux,  the  constancy  of  electron 
spectrum  with  time  and  space,  and  the  brood  wing  of  electron  flux  extending 
hundreds  of  milometers  north  of  the  intense  central  band,  suggest  that  the 
phenomena  leading  to  electron  injection  and  decay  may  be  vastly  different  from 
the  high  altitude  injection  by  trapped  fission  fragments  proposed  earlier. 

The  broad  wing  to  the  north  of  the  band  is  possibly  produced  by  the  decay  of 
fission  fragments  which  were  either  Initially  ulonized  or  which  became  un¬ 
charged  by  charge  exchange  collisions  with  the  atmospheric  atoms.  The  uncharged 
fragments  could  move  freely  across  magnetic  field  lines.  The  existence  of  this 
wing  only  4l  minutes  after  the  detonation  and  the  constancy  of  spectrum  as  a 
function  of  latitude  suggest  that  electrons  could  not  have  diffused  to  this 
region  by  scattering  with  the  atmosphere.  An  alternate  explanation  might  be 
that  the  wing  electrons  were  injected  by  debris  carried  aloft  in  the  rising 
fireball.  However,  in  order  to  place  electrons  in  the  observed  positions,  the 
debris  would  have  to  rise  vertically  almost  1000  km.  Since  debris  rise  could 
produce  a  wing  only  on  the  high  latitude  side  of  the  central  band,  the  occurrence 
of  a  similar  wing  on  the  low  latitude  side  would  Indicate  that  debris  rise  was 
not  the  cause  of  the  wing.  The  data  obtained  with  the  satellite  Explorer  IV 
should  be  examined  closely  for  evidence  of  such  low-latitude  wing. 
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The  time  variation  of  the  electron  flux  is  not  understood,  and  the 
theory  of  electrons  trapped  in  the  geomagnetic  field  is  not  adequate  at 
present  to  give  quantitative  values  of  the  flux  variation  with  time.  Early 
attempts  to  estimate  time  decay  rates  led,  with  certain  extreme  approximations, 
to  a  time  dependence  given  by  l/t.  Another  approach  leads  to  the  conclusion 
that  the  decay  time  is  not  a  simple  function  of  time  but  depends  upon  both 
the  initial  energy  distribution  and  pitch  angle  distribution  of  the  electrons. 
In  essence,  the  distribution  function  of  the  electrons  can  be  decomposed  into 
a  series  of  spatial  eigenfunctions,  each  of  which  decay  in  a  prescribed  manner. 
Only  in  the  simple  case  when  a  single  eigenfunction  is  present  will  the  decay 
rate  be  the  same  at  all  altitudes.  However,  when  the  initial  transients  have 
died  out  and  the  flux  does  decay  with  the  time  constant  appropriate  to  the 
longest  lived  or  dominant  eigenfunction,  the  decay  time  expected  is  on  the 
order  of  tens  of  years,  far  longer  than  the  experimental  results  indicate. 

However,  the  interpretation  of  the  Jason  experiment  in  terms  of  this 
theory  is  difficult  due  primarily  to  the  departure  of  the  geomagnetic  field 
from  that  of  a  centered  dipole.  As  mentioned  previously,  the  trajectories  of 
the  Jason  rockets  lie  in  the  "magnetic  shadow"  of  the  earth.  In  effect,  this 
means  either  that  the  fission  fragments  emitting  the  electrons  must  have  mi¬ 
grated  thousands  of  miles  in  longitude  or  that  the  electrons  seen  by  the  Jason 
vehicles  must  have  been  scattered  sufficiently  to  lower  their  mirror  points 
hundreds  of  kilometers  during  a  single  drift  period  about  the  earth.  No  quan¬ 
titative  calculations  of  the  time  required  for  mirror  point  diffusion  are 
available  to  check  this  latter  possibility,  although  qualitative  arguments 
indicate  that  it  is  possible  providing  the  initial  injection  altitude  is  low 
enough. 

The  observed  decay  rate  in  the  wings  is  sufficiently  close  to  the  time 
variation  of  the  gross  beta  activity  of  fission  fragments  that  one  is  tempted 
to  explain  the  observed  decay  by  assuming  that  the  flux  is  in  some  form  of 
equilibrium  with  the  source  and  that  the  observed  time  dependence  results  from 
the  decaying  activity  of  the  fission  debris  rather  than  from  the  lifetime  of 
individual  electrons.  This  situation  would  occur  if  the  fragments  remained 
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in  a  position  where  they  could  inject  electrons  and  the  electron  mirror 
points  were  sufficiently  low  that  most  of  the  electrons  are  lost  after  one 
drift  period  about  the  earth. 

In  addition  to  explaining  the  time  dependence  of  the  electron  flux,  such 
an  injection  phenomenon  would  also  account  for  the  constant  spectrum  over  the 
time  period  of  the  observation.  However,  the  spectrum  observed  on  flight 
2019  shortly  after  the  detonation  is  still  difficult  to  explain  on  this  basis. 

This  continuous  injection  hypothesis  has  two  fundamental  difficulties. 

The  first  difficulty  is  the  low  injection  efficiency  which  would  probably 
occur  if  the  electrons  were  emitted  at  low  altitudes.  The  altitude  of  the 
bomb  detonation  makes  it  very  unlikely  that  any  of  the  fission  fragments  can 
be  magnetically  trapped  high  above  the  atmosphere,  and  it  is  expected  that 
they  will  quickly  loose  energy  and  come  to  thermal  equilibrium  with  the  atmos¬ 
phere  somewhere  between  100  km  and  perhaps  500  km  depending  on  initial  pitch 
angles.  Since  the  initial  mirroring  point  of  an  electron  cannot  be  higher 
than  the  injection  point,  many  of  these  electrons  will  not  live  long  enough 
to  drift  once  around  the  earth.  Some  of  the  electrons  will  have  their  mirror 
points  raised  by  scattering,  but  until  quantitative  calculations  are  available, 
it  is  difficult  to  estimate  how  effective  this  scattering  might  be. 

The  other  major  difficulty  is  the  requirement  that  the  fission  debris 
be  stationary,  at  least  in  latitude,  for  a  duration  of  several  days.  Although 
data  are  meager,  wind  velocities  above  one  hundred  kilometers  are  expected  to 
be  severed,  hundred  kilometers  per  hour,  a  velocity  which  would  quickly  dis¬ 
perse  the  radioactive  fragments  if  they  were  neutral  and  behaved  as  the  neutral 
constituents  of  the  ionosphere.  If  the  fragments  are  charged,  they  might  be 
held  in  place  by  the  magnetic  field,  at  least  at  altitudes  above  the  exosphere 
where  collisions  with  the  neutral  particles  are  infrequent.  Actually,  a  contin¬ 
ual  change  in  the  degree  of  ionization  of  the  fission  fragments  is  expected: 
During  the  slowing-down  process,  a  tendency  towards  neutralization  by  charge 
exchange  collisions  prevails;  but  after  thermal  equilibrium  with  the  atmospheric 
constituents  is  reached,  fission  fragments  can  be  charged  by  ionizing  collisions 
and  by  the  action  of  the  solar  ultraviolet.  In  fact,  the  fractional  ionization 
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of  the  fission  fragments  would  probably  be  greater  than  that  of  the  atmospheric 
gases  because  of  the  lower  ionization  potentials  of  the  fragment  atoms  com¬ 
pared  to  the  ionization  potentials  of  oxygen  and  nitrogen. 

At  the  present  time,  work  is  continuing  on  these  two  aspects  of  the 
Argus  experiment:  l)  the  ultimate  distribution  of  the  fission  fragments, 
and  2)  the  efficiency  of  low  altitude  injection  as  modified  by  atmospheric 
scattering.  Equations  describing  the  time  behavior  of  the  distribution  function 
of  electrons  trapped  in  the  geomagnetic  field  have  been  derived  and  are  now 
being  programmed  for  solution  by  a  digital  computer.  The  results  of  this  cal¬ 
culation  should  at  least  give  a  qualitative  indication  whether  fragments  sus¬ 
pended  at  a  few-hundred  kilometers  altitude  can  inject  enough  trapped  electrons 
to  agree  with  the  measured  fluxes.  With  regard  to  the  possible  suspension  of 
fission  fragments  at  high  altitude,  it  is  hoped  that  forthcoming  geophysical 
data  may  shed  some  light  on  this  question.  Preliminary  reports  have  indicated 
that  a  high  concentration  of  ionized  atoms  of  meteoretic  material  have  been 
detected  at  high  altitude,  suggesting  that  heavy  ions  may  remain  charged  and 
suspended  for  long  periods  of  time.  Further  evidence  may  be  obtained  from 
studies  of  sodium  ions  detected  in  studies  of  the  night  alrglow.  It  is  hoped 
that  as  a  result  of  this  further  work  it  will  be  possible  to  conclude  whether 
the  short-time  injection  mechanism  or  the  continuous  injection  hypothesis 
best  describes  the  trapping  phenomena  observed  in  the  Argus  experiment. 
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GEOMETRIC  FUNCTIONS  of  DETECTOR  CHANNELS 
(Subscript  Refers  to  Channel  No.) 
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'PAULK  A-l  ( rout'd ) 


GEOMETRIC  FUNCTIONS  of  DETECTOR  CHANNELS 
(Subscript  Refers  to  Channel  No.) 
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DISTRIBUTION 

HEADQUARTERS  USAF 


1  Hq  USAF  (AFCOA),  Wash  25,  DC 

1  Hq  USAF  (AFDAP),  Wash  25,  DC 

1  Hq  USAF  (AFORQ),  Wash  25.  DC 

1  Hq  USAF  (AFRDT),  Wash  25,  DC 

1  Hq  USAF  (AFRDT-NU-1),  Wash  25,  DC 

1  Hq  USAF  (AFCIN),  Wash  25,  DC 

1  Hq  USAF  (AFTAC),  Wash  25,  DC 

1  USAF  Dep  IG  for  Insp  (AFCDI-B-3),  Norton  AFB,  Calif 

1  USAF  Dep  IG  for  Safety  (AFINS),  Kirtland  AFB,  N  Mex 

1  FTD  (Library),  Wright-Patterson  AFB,  Ohio 

1  AFOAR,  ATTN:  RRN,  Major  Munyon,  Bldg  T-D,  Wash  25,  DC 

AFCRL,  L.  G.  Hanscom  Fid,  Bedford,  Mass 

1  ATTN:  Tech  Library 

1  ATTN:  CRZCT,  H.  P.  Gauvin 

1  ATTN:  Dr.  P.  Newman 

1  AFOSR,  ATTN:  SREC,  Bldg  T-D,  Wash  25,  DC 

1  Aeronautical  Research  Lab,  ATTN:  RRLO,  Mr.  Cady,  Wright- 

Patterson  AFB,  Ohio 

MAJOR  AIR  COMMANDS 

AFSC,  Andrews  AFB,  Wash  25,  DC 
I  ATTN:  SCR 

1  ATTN;  SCR -2,  Maj  Berrier 

I  AUL,  Maxwell  AFB,  Ala 

1  USAFIT  (USAF  Institute  of  Technology),  Wright-Patterson  AFB, 

Ohio 

AFSC  ORGANIZATIONS 

1  ASD  (ASAPRL,  Tech  Document  Library),  Wright-Patterson  AFB, 

Ohio 

1  DCAS,  ATTN:  Technical  Data  Center,  AF  Unit  Post  Office,  Los 

Angeles  45,  Calif 

BSD  (Ballistic  Systems  Division),  AF  Unit  Post  Office,  Los 
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Angeles  45, 

Calif 

1 

ATTN: 

BSR-VP-3 

1 

ATTN: 

BSR 

1 

SSD  (Space  Systems  Divison),  ATTN: 
Post  Office,  Los  Angeles  45,  Calif 

ESD,  Hanscom  Fid,  Bedford,  Mass 

1 

ATTN: 

ESAT 

1 

ATTN: 

ESRH 

1 

ATTN: 

ESS 

1 

ATTN: 

ESF 

AFMsl  Dev  Cen  (MDGRT),  Holloman  AFB,  NMex 
AFFTC  (FTFT),  Edwards  AFB,  Calif 

AFMTC,  ATTN:  MU-135,  Tech  Library,  Patrick  AFB,  Fla 
AEDC  (AEOI),  Arnold  Air  Force  Station,  Tenn 
APGC  (PGAPI),  Eglin  AFB,  Fla 
RADC,  Griffis s  AFB,  NY 


KIRTLAND  AFB  ORGANIZATIONS 


AFSWJ,  Kirtland  AFB,  N  Mex 
1  (SWNH) 

48  (SWOI) 

1  (SWT) 

2  (SWRPT,  Capt  Welch) 

1  (SWRPA) 

1  (SWRPI,  Lt  Holly) 

1  (SWRPS,  Maj  Nadler) 

1  (SWRJ) 

1  (SWRB,  Lt  Tinney) 

1  (SWRP,  Maj  Allen) 

1  US  Naval  Weapons  Evaluation  Facility  (NWEF)  (Code  404), 

Kirtland  AFB,  N  Mex 

OTHER  AIR  FORCE  AGENCIES 

l  Director,  USAF  Project  RAND,  ATTN;  RAND  Library,  via:  Air 

Force  Liaison  Office,  The  RAND  Corporation,  1700  Main  Street, 
Santa  Monica,  Calif 
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Mfli  Cyg 

1  Avn  School  of  Medicine,  USAF  Aerospace  Medical  Center,  (ATC), 

Brooks  AFB,  Tex 

ARMY  ACTIVITIES 

1  Chief  of  Research  and  Development,  Department  of  the  Army, 

ATTN:  Lt  Col  Conarty,  Wash  25,  DC 

1  Commanding  Officer,  Diamond  Ordnance  Fuze  Laboratories. 

(ORDTL  06.  33),  Technical  l  ibrary,  Wash  25,  DC 

1  Commanding  Officer,  Chemical  Warfare  Laboratories,  (Library), 

Army  Chemical  Center,  Md 

1  Commander,  Army  Rocket  and  Guided  Missile  Agency  (ARGMA), 

ATTN:  Tech  Library,  Redstone  Arsenal,  Ala 

Director,  Ballistic  Research  Laboratories,  Aberdeen  Proving 
Ground,  Md 

1  ATTN:  Mr.  E.  O.  Baicy 

1  ATTN:  Library 

1  Director,  Evans  Signal  Laboratory,  Weapons  Effects  Section, 

Belmar,  NJ 

1  Commanding  Officer,  US  Army  Signal  Research  &  Development 

Laboratory,  ATTN;  Dr.  W.  McAfee,  Fort  Monmouth,  NJ 

1  Commanding  General,  White  Sands  Missile  Range  (Technical 

Library),  White  Sands,  N  Mex 

NAVY  ACTIVITIES 

1  Chief  of  Naval  Operations,  (OP-36),  Department  of  the  Navy, 

Wash  25,  DC 

1  Commanding  Officer,  Naval  Research  Laboratory,  ATTN:  Tech 

Library,  Wash  25,  DC 

1  Commanding  Officer,  Naval  Radiological  Defense  Laboratory 

(Technical  Info  Div),  San  Francisco  24,  Calif 

1  Office  of  Naval  Research,  ATTN:  Dr.  W.  J.  Thaler,  Wash  25,  DC 

OTHER  DOD  ACTIVITIES 

Chief,  Defense  Atomic  Support  Agency,  Wash  25,  DC 
1  ATTN:  Document  Library 

1  ATTN:  Major  MacCormac 

1  ATTN:  Cmdr  Ennis 

1  Commander,  Field  Command,  Defense  Atomic  Support  Agency 

(FCAG3,  Special  Weapons  Publication  Distribution),  Sandia  Base, 

N  Mex 


LOCKHEED  AIRCRAFT  CORPORATION 


130 


MISSILES  and  SPACE  DIVISION 


TR-61-82 


I 

I 
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Director,  Weapon  Systems  Evaluation  Group,  Room  2E1006, 

The  Pentagon,  Wash  25,  DC 

Director,  Advanced  Research  Projects  Agency,  Department  of 
Defense,  ATTN:  Lt  Col  Werdler,  The  Pentagon,  Wash  25,  DC 

Director,  Defense  Research  fc  Engineering,  The  Pentagon, 

Wash  25,  DC 
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